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PREFACE 

The  David  W.  Taylor  Lectures  were  conceived  to  honor  our  founder  in 
recognition  of  his  many  contributions  to  naval  architecture  and  naval  hydro¬ 
dynamics.  Admiral  Taylor  was  a  pioneer  in  the  use  of  hydrodynamic  theory 
and  mathematics  for  the  solution  of  naval  problems.  He  established  a  tra¬ 
dition  of  applied  scientific  research  at  the  "Model  Basin"  which  has  been 
carefully  nurtured  through  the  decades  and  which  we  treasure  and  maintain 
today.  It  is  in  this  spirit  that  we  have  invited  Professor  Hajime  Maruo 
to  be  a  David  W.  Taylor  Lecturer. 

Prof.  Maruo  was  born  in  Yokohama  in  1922.  He  received  his  profession¬ 
al  education  at  the  University  of  Tokyo,  obtaining  his  Ph.D  in  Naval 
Arc-h If ecture  in  1946.  For  many  years  Prof.  Maruo  has  been  at  the  Yokohama 
National  University  where  he  is  Professor  in  the  Faculty  of  Engineering. 

He  is  a  leader  of  ship  research  in  Japan,  is  an  officer  of  the  Japanese 
Society  of  Naval  Architects,  and  has  been  a  member  of  the  Resistance 
Committee  of  the  International  Towing  Tank  Conference  as  well  as  other 
important  professional  groups.  He  is  an  international  authority  on  ship 
hydrodynamics  and  is  well  known  for  his  very  significant  theoretical  and 
mathematical  research  on  ship  wave  resistance  and  ship  motions  problems. 

Prof.  Maruo  is  no  stranger  to  the  West.  He  spent  a  year  at  Cambridge 
University  in  England  doing  research  under  G.I.  Taylor,  and  in  1964  was  a 
Visiting  Scientist  and  Professor  at  Stevens  Institute  of  Technology.  Over 
the  years  he  has  established  and  maintained  close  ties  with  many  ship 
researchers  in  both  the  United  States  and  Europe. 


FOREWORD 

It  Is  a  great  honor  for  me  to  deliver  the  Sixth  David  W.  Taylor 
Lecture  series  and  I  would  like  to  express  my  sincere  thanks  to  the  David 
Taylor  Naval  Ship  R&D  Center  for  the  kind  arrangement  which  has  enabled  me 
to  have  such  an  opportunity. 

The  subject  which  I  am  going  to  talk  about  is  the  problem  of 
engineering.  It  is  not  the  problem  of  mathematics,  although  recent 
progress  in  technology  depends  greatly  upon  mathematical  theories,  and 
mathematics  has  become  an  indispensable  aid  to  solve  engineering  problems 
nowadays . 


CHALLENGE  TO  BETTER  AGREEMENT  BETWEEN  THEORETICAL  COMPUTATIONS 
AND  MEASUREMENTS  IN  SHIP  HYDRODYNAMICS 

INTRODUCTION 

There  is  a  great  difference  between  the  idea  of  engineering  and 
that  of  mathematics.  The  substantial  importance  in  engineering  is  the 
practical  utility.  Any  theories  cannot  become  useful  unless  they  lead  to 
results  which  are  faithful  representations  of  actual  phenomena. 

In  natural  sciences  like  physics  or  chemistry,  one  may  be  contented 
with  qualitative  agreement  between  theoretical  predictions  and  observa¬ 
tions,  but  in  engineering,  merely  qualitative  agreement  is  not  a  sufficient 
condition  of  the  practical  usefulness.  The  agreement  should  be  quantita¬ 
tive  within  required  accuracy.  Mathematics,  on  the  other  hand,  emphasizes 
logical  rationality.  If  one  considers  an  approximation,  the  rigorous 
mathematical  idea  requires  that  the  simplification  should  be  consistent 
within  itself  throughout  the  approximation.  However,  we  often  encounter 
cases  that  mathematical  rationalism  contradicts  the  practical  usefulness. 

We  know  quite  a  few  cases  that  a  mathematical  theory,  with  all  its  rational 
construction,  yields  rather  unrealistic  results  when  compared  with  the 
actual  phenomena.  On  the  other  hand,  there  are  a  number  of  cases  that 
deviation  from  the  rational  formulation  can  result  in  much  better  agree¬ 
ment  with  measurements.  Generally  speaking,  such  inconsistent  approaches 
are  not  necessarily  safe,  because  their  justification  is  hardly  obtained 
from  the  purely  theoretical  point  of  view.  However,  the  utility  of  this 
kind  of  approximate  method  can  be  appreciated  in  the  practical  appli¬ 
cation,  since  any  theory  which  has  failed  to  give  correct  predictions  is 
almost  useless  even  though  it  has  a  complete  logical  construction  from  the 
mathematical  point  of  view. 

The  present  lecture  intends  to  illustrate  how  the  deviation  from  the 
mathematical  rationality  can  improve  the  agreement  with  measured  results 
and  how  the  mathematical  theory  may  be  revised  for  practical  usefulness. 
Whether  any  inconsistent  approach  can  become  really  useful  or  not  depends 
greatly  upon  the  engineer's  intuition. 

The  title  of  the  lecture  can  cover  very  wide  aspects,  but  I  will 
confine  topics  only  in  problems  of  free  surface  flow,  because  it  is  the 


problem  which  seems  to  be  most  peculiar  to  the  hydrodynamics  of  ships. 

The  lecture  will  deal  with  four  topics,  namely: 

1.  Wave  resistance  of  ships  in  uniform  forward  motion, 

2.  Hydrodynamic  forces  on  oscillating  slender  ships, 

3.  Wave  pressure  on  slender  ships,  and 

4.  Added  resistance  of  ships  in  ambient  ocean  waves. 

WAVE  RESISTANCE  OF  SHIPS  IN  UNIFORM 
FORWARD  MOTION 

Short  Comments  on  Thin  Ship  Theory 

The  theory  of  wave  resistance  is  a  rather  classical  problem.  It  was 

1* 

as  early  as  the  end  of  the  last  century  that  Michell  established  the 

theory  for  thin  ships.  His  theory  was  already  so  complete  that  nothing 

needed  to  be  added  for  a  first  approximation  of  thin  ships.  More  than 

eighty  years  have  passed  since  then;  nevertheless  the  progress  in  the 

theory  of  wave  resistance  has  been  comparatively  slow. 

Michell' s  thin  ship  theory  is  based  on  the  assumption  that  the  beam- 

to-length  ratio  is  so  small  that  its  square  can  be  neglected.  There  are 

several  examples  of  comparison  of  computed  wave  resistance  with  measured 

resistance.  Among  them  are  instructive  results  with  a  series  of  models 

2 

whose  breadth  is  varied  systematically.  It  is  indicated  that  the  wave 
resistance  calculated  by  Michell 's  formula  agrees  well  with  measured 
results,  provided  that  the  beam-to-length  ratio  is  not  greather  than  one 
fifteenth,  as  shown  in  Figure  1.  This  criterion  of  beam-to-length  ratio  is 
too  small  for  practical  hull  forms.  The  beam-to-length  ratio  of  practical 
hulls  is  at  least  one  seventh.  Consequently,  considerable  discrepancy 
appears  between  theory  and  experiment  in  conventional  hull  forms.  This 
fact  does  not  mean,  however,  that  the  thin  ship  theory  is  useless  for 
practical  purposes.  It  is  known  that  the  Michell  thin  ship  theory  has  be¬ 
come  a  very  powerful  tool  for  the  purpose  of  designing  low  resistance  hull 

forms.  A  direct  application  of  the  theory  for  this  purpose  is  known  as 

3 

the  theory  of  minimum  wave  resistance.  The  method  is  an  application  of 
calculus  of  variations.  One  can  determine  the  optimum  curve  of  a  sectional 


*A  complete  listing  of  references  is  given  on  page  119. 
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Figure  1  -  Comparison  of  Michell  Resistance  with  Measured  Wave- 
Pattern  Resistance  (Tsutsumi) 

area  by  which  the  wave  resistance  is  made  minimum  under  appropriate  side 
conditions.  One  very  interesting  fact  is  that  the  curve  of  the  sectional 
area,  which  is  determined  by  the  thin  ship  theory  in  such  a  way  that  the 
wave  resistance  becomes  minimum  at  Froude  number  0.25  with  the  condition 
of  prismatic  coefficient  being  0.60,  is  nearly  identical  with  the  curve  of 

4 

sectional  area  of  the  corresponding  Taylor  Standard  Series  (Figure  2) .  It 
is  recognized  that  the  latter  form  shows  excellent  characteristics  at 
medium  Froude  numbers,  and  this  fact  is  enough  to  testify  the  great  genius 
and  intuition  of  the  late  Admiral  Taylor. 

For  the  purpose  of  quantitative  prediction  of  wave  resistance,  on  the 
other  hand,  Michell' s  theory  is  far  from  useful.  Two  reasons  are  con¬ 
sidered  for  the  discrepancy  of  the  theory  from  measured  results;  one  is  the 
effect  of  viscosity  and  the  other  is  the  effect  of  finite  breadth  of  the 
ship.  The  effect  of  viscosity  upon  wave  resistance  has  been  studied  by 
several  authors,"*  but  no  reliable  theory  has  been  developed  yet.  As  to  the 
effect  of  finite  breadth,  on  the  other  hand,  several  attempts  are  known 


A 
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Figure  2  -  Curve  of  Sectional  Area  of  Minimum 
Wave- Resistance 


which  intend  to  find  out  second  order  terms  with  respect  to  the  beam-to- 
length  ratio  by  means  of  the  successive  method  starting  with  Michell's 
first  order  solution.^  However,  none  of  these  approaches  have  been  proven 
successful  so  the  second  order  theory  is  not  likely  to  be  promising.  We 
can  imagine  the  following  fact  as  the  reason  for  this  difficulty.  It  is 
known  that  the  first  order  solution  is  not  a  uniformly  valid  approximation 
and  higher  order  solutions  show  more  singular  behavior.  It  is  suspected 
that  the  resulting  power  series  with  respect  to  the  beam-to-length  ratio 
which  is  assumed  at  the  beginning  of  the  perturbation  procedure  may  not 
converge  when  infinite  terms  are  taken.  That  means  the  perturbation  ex¬ 
pansion  gives  an  asymptotic  series  without  finite  radius  of  convergence. 

If  it  is  so,  the  addition  of  second  order  terms  does  not  guarantee  any 
improvement  in  accuracy  of  the  approximation. 

Furthermore,  the  existence  of  the  continuous  solution  of  exact 
boundary  value  problem  for  the  fluid  motion  around  the  hull  is  even 


doubtful.  Observations  of  the  white  plume  of  the  wave  crest  around  the 
ship,  may  be  an  indication  of  the  nonexistence  of  the  exact  solution. 

It  is  of  some  interest  to  observe  the  numerical  result  by  Guilloton's 

7  8 

method,  which  shows  a  plausible  improvement  in  agreement  with  measured 

wave  profile  and  wave  resistance.  Unlike  the  consistent  successive 

approximation,  the  various  conditions  involved  are  not  satisfied  at  the 

same  order  of  magnitude.  This  is  because  the  field  equation  is  satisfied 

only  in  the  first  order  while  boundary  conditions  are  partly  satisfied  up 

to  the  second  order.  In  spite  of  such  inconsistencies,  Guilloton's 

method  may  be  regarded  more  useful  than  the  consistent  second  order  theory 

from  an  engineering  point  of  view. 

An  Expression  for  the  Solution  of 
Exact  Nonlinear  Boundary  Value 
Problem 

A  special  feature  of  the  free  surface  flow  is  the  nonlinearity  of 
boundary  conditions.  The  direct  nonlinear  analysis  is  applicable  only  to 
the  simplest  case  such  as  monodirectional  waves  in  a  channel.  Because  of 
its  complexity,  the  only  possibility  of  an  analytical  method  describing 
fluid  motion  around  the  ship  hull  is  the  perturbation  analysis.  As 
mentioned  in  the  preceding  section,  only  the  first  order  approximation  is 
useful  for  practical  cases  because  a  higher  order  approximation  cannot 
improve  the  result  in  most  cases.  The  first  order  solution  is  usually 
obtained  by  the  linearization  of  the  free  surface  condition  at  the  be¬ 
ginning.  However,  there  is  a  possibility  of  giving  a  formal  expression 
for  the  solution  which  satisfies  the  exact  nonlinear  boundary  condition  at 
the  free  surface.  This  expression  is  of  no  use  by  itself  for  the  pre¬ 
diction  of  wave  resistance,  but  it  facilitates  general  discussions  of 
perturbation  analysis. 

First,  we  assume  an  inviscid  incompressible  fluid  with  irrotational 
flow.  A  Cartesian  coordinate  system  with  z-axis  vertically  upward  and  the 
axes  of  x  and  y  is  applied  to  the  fluid  on  the  still  water  surface.  It  is 
convenient  to  employ  dimensionless  quantities.  The  length  scales,  e.g., 
x,  y,  and  z,  are  normalized  by  a  characteristic  length,  £,  which  may  be 


taken  as  the  half  length  of  the  ship,  L/2.  If  the  fluid  velocity  is 

normalized  by  the  ship  speed  U,  then  the  acceleration  of  gravity  is  ex- 

2 

pressed  by  the  dynamic  coefficient  of  gravity  yn  =  g£/U  ,  or  as  a  function 
of  the  Froude  number  y^  =  l/(2Fn  ).  Instead  of  a  ship  moving  in  still 
water,  we  assume  a  uniform  flow  of  velocity  U  in  the  direction  of  x  which 
is  opposite  to  the  motion  of  the  ship.  When  a  ship  hull  is  introduced  at 
a  fixed  position  in  the  flow,  the  flow  field  is  specified  by  a  velocity 
potential  x  +  <J>.  The  disturbance  potential  <j>  is  harmonic  in  the  space 
occupied  by  the  fluid.  Then  the  Laplace  equation 

V2  <t>  =  0  (1) 

is  valid  outside  the  hull  surface  and  below  the  free  surface,  which  is 
expressed  by  the  equation 


z  =  C 


(2) 


The  fluid  velocities  are 


u  =  1  +  9<}>/9x,  v  =  94>/9y,  w  =  9<}>/9z 


(3) 


The  fluid  boundary  is  composed  of  the  wetted  hull  surface,  the  sea  bottom, 
and  the  free  surface.  Since  the  fluid  is  assumed  nonviscous,  the  boundary 
condition  on  the  hull  surface  is  that  the  fluid  velocity  is  tangential  to 
the  hull  surface.  It  is  expressed  by 


3(j>  _  _9x 
3n  9n 


(4) 


where  n  is  the  unit  outward  normal  vector  into  the  hull  surface.  Since 
the  free  surface  is  expressed  by  the  equation  z  =  £,  where  t,  is  a  function 
of  x  and  y,  the  kinematical  condition  on  it  is 

uU  +  vl7'w=0  (5) 
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There  is  a  condition  of  constant  pressure  which  is  expressed  by  the 
Bernoulli  equation. 


j  (u2+v2+w2-l)  +  yQ  £  =  0  (6) 

Then  the  free  surface  elevation  is  given  by 

„  1  2  2  2. 

C  =  -r—  (1-u  -v  -w  )  (7) 

'0 

Keeping  in  mind  the  fact  that  the  boundary  conditions  are  satisfied  at  the 
curved  surface  z  =  5  and  that  £  is  a  function  of  x  and  y,  the  explicit  £ 
in  Equations  (5)  and  (6)  is  eliminated  by  the  substitution  of  partial 
derivatives  of  Equation  (7)  with  respect  to  x  and  y  in  Equation  (5).  The 
result  is 


l(U|^+Vfe  +  W'l^)  (u2’w2+w2>  +  Y0  w  =  0  <8> 

If  we  substitute  the  velocities  by  the  expression  of  Equation  (3),  the 
above  equation  becomes 

l(2|_+v*.v)  l?(,|2+^|  +  Vof  .  0  (» 

This  relation  holds  on  the  unknown  surface  z  =  £,  so  that  the  boundary 
condition  is  quite  nonlinear.  The  usual  way  of  solution  is  the  perturba¬ 
tion  expansion  of  the  boundary  condition  assuming  a  small  parameter  which 
relates  to  the  shape  of  the  hull  surface.  The  first  term  of  the  expansion 
is  the  linearized  solution.  In  order  to  make  such  a  linearized  solution 
valid,  a  condition  of  small  disturbance  is  necessary  which  imposes  a 
restriction  on  the  hull  shape.  Instead  of  the  application  of  the  perturba¬ 
tion  expansion  to  the  boundary  condition,  let  us  seek  a  general  expression 


7 


for  the  solution  of  the  Laplace  equation  with  the  boundary  conditions  in¬ 
cluding  nonlinear  terms.  In  order  to  find  such  an  expression,  we  apply 
Green's  theorem  to  the  velocity  potential  <J>  and  an  appropriate  Green's 
function  in  the  space  bounded  by  a  closed  surface.  Figure  3  shows  this 


Figure  3  -  Coordinate  System  and  Control  Surfaces 


enclosed  surface  to  encompass  the  portion  of  the  hull  surface  in  the  lower 
half  space  S,  a  large  vertical  cylinder  £  surrounding  the  ship,  the  portion 
of  the  horizontal  plane  z  =  0  between  S  and  £,  £q  and  the  portion  of  the 
sea  bottom  inside  the  cylinder  E  .  Then  we  have 


where  G(P,Q)  is  Green’s  function  having  a  simple  pole  at  the  point 
P(x,y,z)  =  QCx'.y'.z'),  P  is  a  point  inside  the  enclosed  space  defined 
above,  Q  is  a  point  on  the  boundary  surface  and  n^  is  the  normal  to  the 
boundary  surface  drawn  inward  at  point  Q.  This  expression  is  valid  pro¬ 
vided  that  the  analytic  continuation  of  the  velocity  potential  in  the 
region  0  >  z  >  C,  is  possible.  This  is  the  only  assumption  for  the  above 
formulation.  Now  we  define  the  Green  function  in  such  a  way  that  it 
satisfies  the  Laplace  equation  in  the  lower  half  space  except  a  point 
P  =  Q,  and  boundary  conditions 


+y 


3x 


,2  '0  3z ' 


G(P,Q)  =0  at  z'  =  0 


(11) 


9n  «<*•«>  ’  0 


on  E, 


(12) 


It  is  assumed  that  the  radiation  condition 


&im  Sx.  G(P,Q)  =  0  or  £im  ^  G(P,Q)  =  0 


(13) 


X'->oo 


is  satisfied  too.  Since  d<p/dn  =  0  and  9G/9n  =  0,  the  contribution  by  the 
integral  on  Eg  vanishes.  If  the  velocity  potential  <)>  is  assumed  to  fulfill 
the  radiation  condition  at  infinite  distance,  the  integral  on  E  decays  out 
as  the  radius  of  the  cylinder  tends  to  infinity.  Since  Eq  is  a  horizontal 
plane,  the  normal  n  is  directed  vertically  downwards.  Then  the  integral 
on  Eq  can  be  transformed  as 


-  j|  [c<p’0>  ^ 


9<HQ)  _ 


*(Q)^£^ 

9nQ 


d  S„ 


=  JJ  G(P,Q)  -  4>(x ' , y ' , z  1 )  9G^tQ) 


dx'dy ' 


z'  =  0 


9 


■  JI  b‘ 


q>  way* 


9z'  Y0  3x'2 


dx'dy 1 


z'=0 


because  of  the  relation  of  Equation  (11).  Integrating  by  parts  twice,  the 
second  term  with  respect  to  x',  yields 


JJ  [o(F.q)  3,K^;y'.2'>  1 


Y0  3x'2 


dx'dy 1 


z'=0 


U, 


Q)  -«x\y\*') 


3x’ 


dy' 


z’=0 


where  is  the  intersection  of  S  and  If 


(i±n  li'l 

\  3x2  Y°  ) 


=  -  $(x,y) 


(14) 


z=0 


then  the  velocity  potential  can  be  written  as 


4>(P)  = 


Mi 


G(P,Q)  5M2>-  -MQ) 

3"o  5nq  j 


d  S„ 


4/ 


dy 


z’-O 


'  ‘•’’h  1 


<Kx\y')  G(P , Q) 


dx'dy' 


(15) 


z'  =  0 
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Since  the  boundary  condition  on  the  free  surface.  Equation  (9)  can  be 
written  as 


£  (2  h  +7»-v)  i7+i2 

z=C 


dz 


we  can  solve  for  the  function  4>(x,y)  as 


■Hx.y) 


i(2i;+v+-7)lv+!2 


(16) 


The  first  term  on  the  right-hand  side  of  Equation  (15)  defines  the  sources 
and  dipoles  with  their  axes  in  the  direction  normal  to  the  surface  and 
distributed  over  the  hull  surface  S.  The  second  term  corresponds  to  a  line 
distribution  along  the  water  line  of  sources  and  x-directed  dipoles,  and 
the  third  term  means  the  source  distribution  over  the  horizontal  plane  or 
the  still  water  plane. 

Next  let  us  show  that  the  same  potential  can  be  expressed  by  a 
distribution  of  sources  only.  The  velocity  potential  given  by  Equation 
(15)  is  valid  outside  the  hull  surface  S.  Here  let  us  assume  a  fictitious 
velocity  potential  <j>  which  is  valid  inside  the  surface  S  and  satisfies  the 
linearized  boundary  condition  at  z  =  0  as 


2  *  * 
0_  +  Y  ii_.  o 
ax2  0  8z 


(17) 


Consider  a  closed  surface  composed  of  S  and  the  portion  of  the  plane  z  *  0 

*  * 
inside  S  denoted  by  Iq.  Apply  Green  s  theorem,  as  before,  to  4>  and 
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G(P,Q)  in  the  interior  domain  of  this  closed  surface.  If  the  point  P  is 
outside  S,  we  have 


1_ 


f  f  G(P,Q)  ^  ^  -<f>*(Q) 
^S+I*  L  3nQ 


9G(P,Q) 
dnl  J 


d  S, 


0 


where  n  is  the  normal  drawn  inwards  to  the  domain  under  consideration. 
Integrating  by  parts  over  £q,  as  before,  we  find 


Subtraction  of  Equation  (18)  from  Equation  (15)  yields 


<f>(P) 


■  -  fcff  + “ + *‘(Q)  “I 


‘^O  I 


G(P,Q)  ^ 


( 34>(Q) 


i Lm 

9x' 


{<KQ>-V'(Q)} 


dy’ 


2  '  =  0 


-Ajl. 


<Kx',y')  G(P,Q)  dx'dy' 

z'=0 


(r 


Now  we  assume  that  the  fictitious  velocity  potential  is  chosen  in  such  a 
way  that  it  has  a  value  identical  with  <J>  on  S.  Then,  the  relations  on  S 


0 


,  3G  ,  *  3G 


,  3G  *  3G 

=  *  IK  ~  ♦  3^  = 


«jHO  f  =  o 


Therefore,  Equation  (19)  becomes 


d  S, 


4 1  b 


o)  (mm  _  m!m  | 

g'  |  3x'  3x'  j 


dy* 


z  -0 


4tty  II  $<x *y')  G(P»Q)  dx'dy ' 
’"W  lz'=0 

L0 


(20) 


Since  -G(P,Q)  means  a  source  at  the  point  Q,  the  velocity  potential  is  ex¬ 
pressed  by  the  distribution  of  sources  over  the  hull  and  the  still  water 
surface.  The  density  of  the  hull  surface  sources  is 


Here  we  take  a  plane  parallel  to  x  involving  the  normal  and  take  a  length 
s^  along  the  curve  of  intersection  of  this  plane  and  the  surface  S.  If  a 
is  the  angle  between  the  normal  and  the  x  axis,  we  have  the  relation 

3$  3$  3<j> 

-r* -  =  cos  a  +  sin  a  -r-1-— 

9x  3n  3s^ 

* 

We  have  assumed  that  <j>  =  <J>  on  S,  so  that 
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Then  the  line  integral  is  expressed  by 


J  G(P,Q)  ( 


|i_|Ljcos  ady'  =  J 


G(P,Q) 


(fi* 


n  ^-ds 
3n*/  X  ds 


=  4u  J*  o(Q)  G(P,Q)  nx  ds 


where  s  is  the  length  along  and  n^  =  3x/9n.  Then  we  obtain 


4>(P) 


‘  ‘  II 0 


(Q)  G(P,Q)  d  SQ  + 


ki. 


o(Q)  G(P,Q)  nx  ^  ds 


4><x\y')  G(P,Q)|  ^  ^  dx'dy’ 


(21) 


If  iwe  assume  the  disturbance  velocity  is  so  small  that  the  function  $(x,y) 
is  a  negligible  second  order  contribution,  the  velocity  potential  can  be 
expressed  by  the  source  distribution  on  the  hull  surface  accompanied  by 
the  line  distribution.  This  yields 


<KP) 


G(P»Q)  O(Q)  d  SQ 


G(P,Q)  o(Q)  n 


dzl 

ds 


ds 


(22) 


One  can  determine  the  source  density  o(Q)  so  as  to  satisfy  the  boundary 
condition  on  the  hull  surface.  Brard  named  this  kind  of  boundary  value 


problem  the  Neumann-Kelvin  problem.  Several  numerical  works  have  been 
carried  out  so  far. ^  Since  the  smallness  of  the  disturbance  velocity  can¬ 
not  be  assumed  apriori,  other  conditions  are  needed  in  order  to  realize  it. 

The  simplest  case  is  the  thin  ship.  If  the  beam  of  the  ship  is  very  small, 

* 

the  inner  region,  in  which  the  potential  <j>  is  defined,  shrinks  to  a  narrow 

*  * 

slit  and,  consequently,  3<}>  /3n  becomes  higher  order  with  respect  to  the 
beam-to-length  ratio  of  the  ship.  Therefore,  o(Q)  is  determined  by 
(1/ Att)  3t})/3n.  The  line  integral  becomes  third  order  and  can  be  omitted. 

If  we  write  the  equation  of  the  hull  surface  as 

y  =  f(x,z)  sgn  y  (23) 

the  velocity  potential  is  reduced  to 

<!>(P)  =  -  ^  JJ  G(P,Q)  fx(x',z')  dx'dz'  (24) 

S 

c 

where  is  the  center  plane  of  the  ship  and  =  3f/3x.  This  is  identical 
with  Michell's  potential. 

Wave  Resistance  at  Low  Speed 

Most  existing  theories  of  ship  waves  and  wave  resistance  are  based  on 
the  linearization  of  the  flow  field  by  a  small  parameter  which  specifies 
the  slenderness  of  the  ship  hull.  Since  ship  hulls,  in  practice,  are 
neither  so  slender  nor  thin  enough  to  secure  the  validity  of  the  linearized 
theory,  the  agreement  between  the  theoretical  prediction  and  the  experi¬ 
mental  results  is,  in  general,  not  satisfactory.  The  disturbance  ve¬ 
locities  are  not  small  enough  to  make  their  square  negligible  everywhere 
on  the  free  surface.  Since  the  inclusion  of  the  nonlinear  terms  in  the 
free  surface  condition  makes  the  boundary  value  problem  intractable,  some 
simplification  other  than  the  linearization  by  the  beam-to-length  ratio 
as  a  small  parameter  is  needed  to  formulate  the  wave  resistance  of 
practical  hull  forms,  expecially  full-form  ships. 


Because  the  operating  speed  of  ordinary  ships  of  displacement  type 
is  in  the  range  of  Froude  numbers  from  0.15  to  0.30,  one  may  regard  that 
usual  ships  are  operating  in  comparatively  low  velocity.  If  the  speed  of 
advance  is  extremely  low  and  the  elevation  of  the  free  surface  is  very 
small,  the  flow  around  the  hull  is  comparable  to  flow  with  an  undisturbed 
free  surface,  and  is  similar  to  the  flow  around  a  double  body  fixed  in  a 
uniform  stream.  Then  the  deviation  of  actual  flow  from  the  double  body 
flow  is  due  to  the  elevation  of  the  free  surface.  Since  the  elevation  of 
the  free  surface  depends  on  the  Froude  number,  we  may  employ  the  pertur¬ 
bation  expansion  of  the  free  surface  condition  by  the  Froude  number  as  a 
small  parameter,  but  it  will  be  found  later  that  this  approach  is  not  so 
simple. 

If  our  purpose  is  to  formulate  the  wave  resistance  at  low  Froude 
numbers,  we  can  derive  an  approximate  formula  directly  from  a  general 
expression  of  the  wave  resistance.  The  wave  resistance  is  determined  by 
the  momentum  or  energy  analysis  of  the  asymptotic  expression  of  the  fluid 
motion  at  a  great  distance  from  the  ship. 

The  Green's  function  G(P,Q)  has  an  asymptotic  expression  when  the 
point  P  is  brought  to  infinite  downstream  x  -►  00 .  If  we  assume  the  case  of 
infinite  depth  of  water,  it  takes  the  form 


G(P,Q) 


Yo(z+z')sec^0 

e 


sinCy^x-x'  sec  0) 


cos(Ygy-y'  sec 


0  tan  0)  sec  0d0  (25) 
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The  fluid  motion  there  is  characterized  by  the  Kochin  function.  In  the 

13 

case  of  the  distributed  Havelock  sources  of  density  c(x,y,z)  on  the 
surface  S,  such  as  the  first  term  of  the  right-hand  side  of  Equation  (22), 
the  Kochin  function  becomes 


JI 


H(k,0)  =  - 


o(x,y,z)  exp[kz+ik(x  cos  0  +  y  sin  0)]  dS 


(26) 


In  the  case  of  Equation  (15),  on  the  other  hand,  the  singularities  are  the 
surface  distribution  of  sources  and  dipoles  on  S,  the  line  distribution  of 
sources  and  dipoles  along  L,  and  the  distribution  of  sources  on  the  hori¬ 
zontal  plane  E^.  Then,  the  corresponding  Kochin  function  is 

H(k, 0)  =  -  ~  (-|^  -  )  exp[kz+ik(x  cos  6  +  y  sin  0)  ]  dS 


ik  cos  0 


exp[ik(x  cos  0  +  y  sin  0)]  dy 


4>(x,y)  exp[ik(x  cos  0  +  y  sin  0)]  dxdy 


In  order  to  distinguish  the  above  function  from  the  ordinary  Kochin 
function  such  as  Equation  (26),  it  may  be  called  the  generalized  Kochin 
function.  Keeping  in  mind  the  fact  that  4>(x,y )  decays  out  much  faster  than 
the  fluid  velocity  on  going  away  from  the  ship,  one  can  express  the  wave 
resistance  by  Havelock's  formula 


2  2  2 

R  =  8irPU  A 
w  '  0 


sec^G,  0)|^sec^0  d0 


Now,  let  us  consider  the  first  approximation  for  low  Froude  numbers.  It  is 

easily  understood  by  the  condition  of  Equation  (5)  that  the  vertical 

2 

velocity  w  at  the  free  surface  is  of  the  order  of  Fn  ,  since  the  free 

2 

surface  elevation  is  0(Fn  )  by  Equation  (7).  Therefore,  the  zeroth 
approximation  for  the  velocity  potential  at  low  Froude  numbers,  designated 
by  <{)q,  is  obtained  from  the  condition 


3<|>q/9z  =0  at  z  =  0 
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The  second  and  third  terms  on  the  right-hand  side  have  the  factor  or 

2  ° 

Fn  ,  while  the  first  term  does  not,  so  that  one  may  consider  the  former 

two  terms  are  of  higher  order.  This  is  not  the  case  because/  the  first  and 

second  terms  can  be  transformed  into  an  integral  over  the  plane  z  =  0.  To 

show  this,  let  us  express  the  double  body  potential  in  the  form 


Vp>  '  - 


//  .[ 
s+s 


34>0(Q)  1 

-  ~  -  VQ) 

PQ 


8n„ 


d  S. 


(34) 


where  S  is  the  mirror  image  surface  of  S  in  reference  to  the  plane  z  =  0 

- -  * 

and  PQ  is  the  distance  between  P  and  Q.  If  the  point  P  is  inside  S+S  , 

2  2 

then  <j>g(P)  =  0.  Next  we  consider  the  Fourier  transform  of  3  4>q/3x  on  the 
plane  z  =  0.  Since  4>0  =  0  inside  S  +  S  ,  we  have,  after  integrating  by 
parts  twice  with  respect  to  x 


£ik(x  cos  8  +  y  sin 


6) 


dxdy 


,  2  2 
=  -  k  cos 


II 


ik(x  cos  0  +  y  sin  8) 


dxdy 


-j  e*- 


ik  cos  0*0, 


gik(x  cos  0  +  y  sin  0)  ^ 


(35) 


z=0 


The  value  of  on  the  plane  z 


0,  on  the  other  hand,  becomes 
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1_ 

2tt 


34>0(Q)  ! 

=  -v<» 


3n 


Q  PQ 


-(- 

Q  \PQ 


z=0 


d  S, 


(36) 


because  of  the  symmetry.  Now  we  employ  the  integral  representation 


TT 


—  =  J*  d0  exp[-k{|z-z' |+i(x-x')cos  0  +  i(y-y’)sin  0}]  dk  (37) 


in  Equation  (34).  Then  the  Fourier  transform  of  4>q)  becomes 

z=0 


00  . 

II  V2=0  exp[ik( 


x  cos  0  +  y  sin  0)]  dxdy 


34>0(Q) 


3n„ 


'*o(Q)  ^ 


exp[kz’+ik(x'  cos  0  +  y' sin  0)]  d 


(38) 


Combining  Equations  (35)  and  (38),  we  find 


h  JJ  (dl r  -*0 1^)  e*ptkz+ik<x 


cos  0  +  y  sin  0) ]  ds 


,-L-  f  (*Z 
4"Yo  J.  '  ;)x 


-ik  cos  0*$  )  exp[ik(x  cos  0  +  y  sin  0))  dy 

U/z=0 


1  3  ^0  k  /.  k  cos20 \ ~ 

4’Y0  5*2  4”  V  Yo  / 


(39) 
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Insert  the  above  relation  in  Equation  (32)  and  put  k  *=  Yq  sec  6.  If  we 

write  the  relative  flow  velocity  around  the  double  body  as  un,  vn>  w.,  and 

2  2  2  u  u  u 

put  Uq  +  Vq  *  q^,  after  some  reductions,  we  have  the  expression  as 


H(Y0  sec20,0)  =  "  a^T  JJ  ^  S(q0‘1)}  +  {v0(V1>: 


z=0 


expfiy^  sec  0(x  +  y  tan  0)]  dxdy 


■Ml 


D(x,y)  exp[iyQ  sec  0(x  +  y  tan  0)]  dxdy  (40) 
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where  the  function  D(x,y)  was  defined  first  by  Baba  and  is  written  as 


D(x>y)  +  3y  ^0V0^ 


(41) 


We  may  derive  another  expression  from  the  source  distribution  over  the 
hull  surface  S  given  by  Equation  (21).  Then  the  generalized  Kochin 
function  becomes 


-  II « 


H(k,0)  =  -  |  I  o(x,y,z)  exp(kz  +  ik(x  cos  0  +  y  sin  0)]  dS 

S 


h~  o(x,y, 

Y0  JT 


0)  exp[ik(x  cos  0  +  y  sin  0)]  n^  dy 


■  4”Yo  11 


$(x,y)  exp[ik(x  cos  0  +  y  sin  0)]  dxdy  (42) 
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The  first  approximation  at  low  speed  is  given  by  the  substitution  of 
o(x,y,z)  by  a0(x,y,z)  and  4>(x,y)  by  4>0(x,y),  where  a0(x,y,z)  is  the  source 
distribution  of  the  double  body  in  a  uniform  flow.  The  result  of  this 
approximation  is  not  identical  with  that  of  Equation  (32)^  because  the 
boundary  condition  on  the  hull  surface  is  not  satisfied  by  the  source 
distribution  Og(x,y,z). 

As  a  numerical  example,  wave  resistance  of  Wigley's  parabolic  model  is 
calculated. ^  Figure  4  shows  the  result  calculated  by  Equation  (28)  with 
the  Kochin  function  defined  by  Equations  (32)  or  (40)  and  by  (42).  The 
results  are  compared  with  the  residuary  resistance  of  towing  tests  and 
results  of  wave  pattern  analysis  of  the  longitudinal  cut  method  as  well  as 
the  calculation  by  Michell's  formula.  A  considerable  improvement  is 
observed  in  agreement  with  the  measured  results,  especially  in  lower  speed 


x  10 


0.15 


0  20 


0  25 


u  /  )V 


0  30 


0  35 


Figure  4  -  Computed  and  Measured  Wave-Resistance  Coefficient 

of  Wigley  Model 
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by  the  computation  with  Equation  (32).  However,  there  still  exist  some 
discrepancies  when  the  Froude  number  is  greater  than  0.20. 

In  order  to  achieve  further  improvement,  some  kind  of  revision  of  the 
formula  is  desirable.  We  have  employed  the  double  body  source  distribution 
in  the  expression  of  Equation  (42),  but  it  does  not  satisfy  the  hull  sur¬ 
face  boundary  condition  at  finite  Froude  numbers.  Therefore,  we  try  to 
employ  the  source  distribution  by  which  the  hull  surface  boundary 
condition  is  satisfied  at  the  Froude  number  under  consideration.  If  we 
assume  that  the  line  distribution  along  the  water  line  and  the  distribution 
over  the  horizontal  plane  do  not  make  serious  effects  on  the  hull  surface 
boundary  condition,  the  density  of  the  source  distribution  over  S  is 
determined  by  the  integral  equation 


n  =  -  27ra(P)  + 
x 


a(Q)  G(P,Q) 
p 


d  S„ 


(43) 


A  numerical  method  is  available  to  determine  the  source  density  o(P). 

Then  the  result  is  substituted  in  the  formula.  Equation  (42) ,  to  calculate 
the  generalized  Kochin  function,  which  determines  the  wave  resistance  by 
the  formula.  Equation  (28).  A  numerical  example  is  given  for  the  case  of 
Yq  =  6.0  or  Fn  =  0.2887.  Three  kinds  of  calculations  are  compared  in 
Figure  4.  The  first  is  the  calculation  by  the  ordinary  Kochin  function  of 
the  distribution  of  sources  on  S  only,  i.e.,  only  the  first  term  on  the 
right-hand  side  of  Equation  (42)  is  taken.  The  second  is  the  addition  of 
the  line  integral,  which  corresponds  with  the  result  of  the  Neumann- 
Kelvin  solution  suggested  by  Brard.  The  third  example  is  the  inclusion  of 
all  terras  of  Equation  (42).  Much  closer  agreement  with  measured  results 
are  obtained  by  the  addition  of  the  third  term  of  Equation  (42).  Recently 
Kitazawa^  et  al.,  presented  another  approximation.  They  put  the  density 
of  hull  surface  sources  as 


o(P)  =  oQ(P)  +  o1(P) 

and  determined  o^(P)  by  the  integral  equation 


(44) 
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-  °1(P)  +  S  JJ  °1«>  3^  G(P’Q>  d  SQ  -  2 ¥^J  °1«>  sf  G<P'«  V*' 


'  si2  II, 


D(x,y)  -rr-  G(P,Q)  dx'dy' 


9n 

p 


(45) 


This  is  equivalent  to  the  velocity  potential 


4><P)  = 


II. a 


(Q)  G(P,Q)  d  SQ  + 


k  I.  0 


(Q)  G(P,Q)  nxdy' 


'  ^ 0  II 


$0(x',y')  G(P,Q)  dx'dy' 


(46) 


where  the  hull  surface  source  density  is  so  determined  that  the  hull  sur¬ 
face  boundary  condition  is  satisfied.  Therefore,  the  only  approximation 
is  the  replacement  of  $(x,y)  by  $g(x,y).  Numerical  results  show  a 
plausible  agreement  with  measured  results  as  shown  in  Figure  4.  However, 
none  of  these  results  are  regarded  as  consistent  approximations  from  the 
rigorous  aspect  of  the  perturbation  analysis. 

Method  of  Coordinate  Transformation 

In  the  preceding  section,  we  assumed  the  deviation  of  the  flow 
velocities  around  the  hull  from  those  of  double  body  flow  as  a  small 
perturbation  and  formulated  the  first  approximation  for  the  wave  re¬ 
sistance.  If  we  want  expressions  for  the  flow  velocity  or  wave  pattern, 
however,  the  perturbation  expansion  of  the  velocity  potential  is  needed. 
Now  let  us  write  the  fluid  velocity  around  the  hull  in  the  form 

u  =  u0  +  ui»  v  =  v0  +  vi*  w  =  wq  +  wi  (47) 
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The  first  approximation  for  the  free  surface  elevation  is  given  by 
Equation  (31)  or 


C 


0 


2  2^ 

(1_uo'vo) 


2=0 


(48) 


Then  the  kinematical  condition  of  the  free  surface.  Equation  (8),  gives 


W,  =  u„  -r — 
1  0  dX 
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at  z=0 


(49) 


This  relation  provides  the  boundary  condition  at  the  free  surface  for  the 

first  approximation  of  the  perturbation  velocity  potential.  However,  the 

right-hand  side  is  determined  by  the  double  body  flow  which  does  not 

present  a  wavelike  motion,  so  that  the  boundary  value  problem  gives  the 

solution  which  is  not  wavelike.  This  result  contradicts  the  actual 

phenomena.  In  order  to  avoid  this  contradiction,  one  has  to  revise  the 

basic  assumption.  Here  we  employ  the  hypothesis  which  was  proposed  by 
18 

Ogilvie.  The  basic  assumption  is  that  the  perturbation  velocity  is 

wavelike  and  the  wavelength  is  proportional  to  the  Froude  number  squared. 

This  means  that  the  differentiation  results  in  the  change  of  order  of 

magnitude  by  the  order  of  wave  number.  Secondly,  it  is  assumed  that  the 

wavelike  nature  appears  in  the  first  approximation  of  the  perturbation 

2 

velocities,  i.e.,  u^,  v^,  w^  which  are  0(Fn  ),  but  not  in  the  first 
approximation  of  the  free  surface  elevation  r,^.  The  second  approximation 
for  the  surface  elevation  is  given  by 


’1 


(50) 


4  2 

and  is  0(Fn  ).  Since  u^,  v^,  w^  are  0(Fn  )  and  their  differentiation 

_2 

reduces  the  order  by  Fn  ,  the  derivatives  of  the  perturbation 
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velocities  are  0(1).  Inserting  the  expression,  Equation  (47)  in  the  free 
surface  condition 


1/  3  3  3  \  .  2  2 .  2.  , 

2  (u  H  ^  ^  (U  ^  +W  >  +  Y0  W  =  ° 


(51) 


and  collecting  terms  of  0(1) ,  we  obtain  on  the  surface  z  =  £ 
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(52) 


The  function  D(x,y)  is  expressed  by  Equation  (41).  If  the  double  body  flow 
velocity  is  known,  the  function  D(x,y)  is  a  known  function.  Then  Equation 
(52)  gives  the  boundary  condition  at  the  free  surface  for  the  perturbation 
velocity  potential  4^  defined  by 


U1  =  3x~ 


3$, 


vi  ay 


3<J>j 

wi  =  dT 


(53) 


as  follows. 


2  324>x  a\  2  a241  3^ 

0  9x2  +  2  u0  V0  3x3y  +  v0  9y2  +  Y0  3z“  =  T0  D(x,y) 


(54) 


The  perturbation  velocity  potential  is  harmonic  outside  the  hull  surface 
S  and  below  the  curved  surface  z  =  £  on  which  the  boundary  condition. 
Equation  (54),  is  satisfied.  Since  the  boundary  condition  on  the  curved 


surface  is  not  convenient  to  finding  the  solution,  let  us  perform  the 
following  transformation  to  the  vertical  coordinate. 

z  =  +  z'  (55) 

After  the  transformation  of  the  Laplace  equation  to  new  coordinates  and 
taking  only  the  first  order  terms,  we  find  that  the  Laplace  equation  does 
not  change,  namely 
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a  a 
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(56) 


while  the  boundary  condition  at  the  free  surface  is 


2  a\  2  a2^  a*! 

0  9x2  +  2  u0  V0  3x3y  +  V0  +  Y0  Jz7  =  Y0  D(x,y) 


(57) 


at  z'  =  0.  The  homogeneous  equation,  which  is  obtained  by  setting  the 
right-hand  side  equal  to  zero,  is  the  boundary  condition  of  free  waves  on 
a  nonuniform  flow  field. 


2  3\  2  d(pl 

U0  T~2~  +  2  U0  V0  3x3y  +  V0  “  2  +  Y0  dz7  =  ° 


(58) 
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Now  let  us  assume  the  wave  potential 


(J>1  =  A  exp  Y0(k(x,y)z'+i  S(x,y)] 


(59) 


On  substituting  in  the  Laplace  equation  and  taking  terms  of  the  lowest 
order  with  respect  to  Yq\  we  have  the  relation 
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Because  of  the  exponential  decay  in  the  z  direction,  the  wave  motion  is 

-1  2 

significant  only  in  a  surface  layer  of  thickness  0(Yn  ).  Therefore,  z'  is 

4  0 

0(Fn  )  and  the  second  term  of  the  above  equation  is  to  be  omitted.  Then 

we  have 


- \/(i)2+(f)2 

that  means  the  local  wave  number.  Substituting  Equation  (59)  in  Equation 
(58)  and  taking  terms  of  the  lowest  order,  we  obtain 


/  3S 
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or 
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This  is  identical  with  the  dispersion  relation  which  has  been  indicated  by 
19 

Keller.  The  solution  of  the  above  differential  equation  determines  the 
phase  function  S(x,y).  Keller  has  proposed  a  kinematical  theory  of  waves 
superimposed  on  the  nonuniform  flow  around  the  ship  by  the  analogy  with 
the  geometrical  optics. 

Now  let  us  show  the  possibility  of  giving  an  analytical  expression  for 
the  wave  function  by  a  coordinate  transformation.  To  find  the  solution, 
we  employ  the  curvilinear  coordinates  along  streamlines  and  equipotential 
lines  of  the  double  body  flow  in  the  plane  z  ■  0.  Designate  the  velocity 
potential  of  the  double  body  flow  by  4>  and  the  stream  function  on  the 
plane  z  -  0  by  ?,  the  latter  of  which  is  slightly  different  from  the  stream 
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function  in  the  two-dimensional  motion.  They  are  defined  by  the  relation 
with  the  velocity  of  the  double  body  flow  as  follows. 


=  30  =  1  35 
U0  3x  h  3y 


(64) 


=  3$  =  _  1  3? 

V0  3y  h  3x 


(65) 


If  the  double  body  is  the  body  of  revolution,  then  h  is  the  radial  length. 
In  general  cases,  it  is  determined  by  the  partial  differential  equation 


111  a.  _ 

U0  3x  V0  3y 


dx  dy 


(66) 


If  we  change  the  independent  variables  of  Equation  (58)  from  x,y  to  $,5, 
there  is  a  relation 


3(f)  3<}>i  ?  30. 

u0  3^r  +  v0  3T  =  q0  3T  (67) 

Because  of  the  fact  that  differentiation  of  the  perturbation  potential 

reduces  the  order  of  magnitude,  the  order  of  30/30  and  30  /35  are  higher 
2  2-1  1 

than  that  of  3  0^/30  by  .  Differentiating  the  above  equation  again 
and  omitting  higher  order  terms,  one  obtains  from  Equation  (57) 

4 

qorr  +  YoF  =  Yo  D(x’y)  (68) 

3$ 

Next,  the  independent  variables  are  transformed  again  to  new  variables  f,, 
H,  C  by  the  following  relations. 


dC 


(69) 


29 


h  dn 
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(71) 


where  q^^.ri)  =  qg(x,y)  and  the  lower  limit  of  integrals  depends  on  the 
origin  of  $  and  4\  If  we  assume  q^  to  be  slowly  varying,  so  that  9  qg/9<l> 
and  9  q^S?  are  0(1),  we  can  prove 
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(72) 


Differentiating  again  by  £  and  taking  only  the  term  of  the  lowest  order, 
we  have 


A  similar  assumption  is  applied  to  9<^/9£  such  as 


34>x 

ar 
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Multiplying  qg/Yg 


on  both  sides  of  Equation  (69)  and 


writing 


2 

q0 

—  D(x,y)  =  E(C,n) 


(75) 


we  obtain 


30 


(76) 
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A  similar  discussion  is  applied  to  the  transformation  of  the  Laplace 
equation,  and  the  lowest  order  term  becomes 

a2$  a24>  a2<t> 

— y-  +  ~^r  +  — £  =  0  <77> 

a?  an  a^ 

Thus,  the  Laplace  equation  keeps  its  original  form  after  the  transforma¬ 
tion.  To  find  the  solution  which  satisfies  the  inhomogeneous  boundary 
condition  of  Equation  (76),  we  consider  the  basic  solution  which  satisfies 
the  boundary  condition 


a24>  a# 

^r+ar 


6(C-S\n-n') 


(78) 


where  6(5-?',n-n')  is  the  delta  function  of  two  variables.  The  solution 
which  satisfies  the  above  at  C  =  0  and  the  radiation  condition  at  infi¬ 
nite  distance  as  well,  is  given  by  the  function 


G(C,n,c;C',n') 
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where  I  means  the  imaginary  part  is  taken.  Then  the  solution  of  the 
boundary  condition  of  Equation  (76)  is  given  by 


+! 


e(5',o')  g^.o.c^’.o')  dC'dn' 


(80) 


The  asymptotic  form  far  downstream,  i.e.,  5  -*■  °°,  is  written  in  the  form 

tt/2 

E(C' »Tl,)dC'dr|,I  exp [ C  sec2  0+i  sec  0(5-5'+  tan  0(0-0'))] 

-tt/2 

•  sec2  0  d0  (81) 


Thus,  the  phase  function  of  the  elementary  wave  defined  by  Equation  (59)  is 

S  =  Yq1  sec  0  (5+0  tan  0) 


Now  let  us  show  that  the  phase  function  satisfies  the  dispersion  relation 
given  by  Equation  (63).  There  are  relations 
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Since  it  can  be  shown  that  3qQ/3o  =  0(Yq  ),  we  have,  after  deleting  higher 
order  terms, 
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A  similar  argument  is  applied  to  3£/3y,  Sr)/9x,  and  Sr|/9y>  as 


_  Yovo  an  Yovo  an  Youo 

ay  3  ’  ax  3  ’  3y  “  3 

J  r\  n  *  n 


Then  we  obtain 
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It  is  readily  shown  that  the  dispersion  relation  of  Equation  (63)  is 
satisfied  by  the  substitution  of  the  above  relations. 

The  amplitude  function  of  the  free  wave  at  infinite  downstream  is 
expressed  by  the  complex  form 


E(£,n)  exp[i  sec  0(£+n  tan  0)]  d£dn 


•-Hf 


D(x,y)  exp[i  sec  0(£+n  tan  9)]  dxdy 


The  integral  with  respect  to  x,  y  is  carried  out  over  the  plane  z  *  0 
excluding  the  interior  of  the  hull. 
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The  second  approximation  of  the  free  surface  elevation  is  given  by 


because  qQ  1  at  a  great  distance. 

The  wave  resistance  experienced  by  the  ship  is  derived  from  momentum 
or  energy  analysis  of  free  waves  at  a  great  distance.  It  is  determined  by 
the  amplitude  function 

tf/2 

R  =  p  U2  H2  | A(0) 1 2  sec3  0d0  (85) 
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In  order  to  calculate  A(0) ,  we  need  the  inverse  transform  of  Equations  (69) 
and  (70).  For  this  purpose,  we  take  curvilinear  coordinates  along 
streamlines  and  equipotential  lines  along  which  the  length  s  and  t  are 
taken,  respectively.  It  is  proved  that  we  have  the  relations 


d^  =  >0rf +  0<^1)1 
K  j 


dn  =  Y0j3  +  0(Yo1:)j 


(86) 


d$  =  ds  dY  =  q^  h  d  t 
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Omitting  higher  order  terms  in  Equation  (86)  and  integrating,  we  obtain 


S  =  Y, 


r  ds  r  dt 
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(88) 
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Since  Yq/<1q  corresponds  to  the  local  wave  number,  the  new  coordinates  mean 

the  streamline  coordinates  with  scales  vary  proportionally  to  the  local 

wave  length.  The  above  result  takes  account  of  the  distortion  of  the  wave 

pattern  due  to  the  nonuniform  base  flow  velocities  near  the  hull.  The 

boundary  condition  on  the  hull  surface  is  satisfied  in  the  case  of 

symmetric  flow,  because  the  velocities  u  and  v  are  tangential  to  the  hull 

surface  at  z'  =  0  and  are  significant  only  in  a  thin  layer  of  thickness 

0(yq1)  near  the  free  surface.  Though  the  theory  developed  here  looks  like 

a  reasonable  representation  of  the  actual  phenomena,  no  numerical  results 

have  been  presented  so  far.  It  should  be  noted  that  a  purely  numerical 

20 

method  has  been  employed  in  a  similar  boundary  value  problem  by  Dawson 
and  the  result  shows  a  plausible  agreement  with  experiments. 


HYDRODYNAMIC  FORCES  ON  OSCILLATING 
SLENDER  SHIPS 

Application  of  the  Slender  Body 
Theory 

As  mentioned  in  the  preceding  chapter,  the  possibility  of  mathematical 
analysis  of  the  fluid  motion  around  a  ship  hull  depends  substant ially  on 
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the  linearization  of  boundary  conditions  at  the  free  surface.  In  the  case 
of  ships  or  other  floating  bodies  making  oscillations  on  the  free  surface 
with  its  average  position  at  a  fixed  point,  indicating  no  average  veloc¬ 
ity,  the  boundary  condition  at  the  free  surface  can  be  linearized  in  a 
simple  fashion  by  assuming  that  the  amplitude  of  the  oscillation  is 
sufficiently  small.  Therefore,  no  restriction  is  imposed  on  the  shape  of 
the  body.  By  the  employment  of  the  numerical  method,  one  can  calculate 
hydrodynamic  forces  on  any  kind  of  shapes  of  floating  bodies  in  principle, 
and  it  is  known  that  some  numerical  results  show  fairly  good  agreement  with 
measured  values.  However,  a  rational  development  of  the  linearized 
theory  becomes  too  intricate  when  the  steady  forward  speed  is  introduced 

to  the  oscillating  ships.  The  difficulty  in  finding  a  rational  solution 

21 

which  is  not  trivial  was  first  demonstrated  by  Peters  and  Stoker.  They 

showed  that  the  hydrodynamic  reactions  such  as  the  added  mass  and  damping 

did  not  appear  in  the  order  of  approximation  of  the  linearized  theory  for 

a  thin  ship  oscillating  in  the  plane  of  symmetry.  There  had  been  extensive 
22  23 

works  by  Haskind  and  Hanaoka  about  thin  ships  in  longitudinal  os¬ 
cillations  in  still  water  before  that  time.  A  full  condemnation  of  these 
achievements  by  the  reason  of  inconsistency  may  be  unfair,  because  the 
consistent  structure  of  theory  breaks  down  on  account  of  just  a  single 
reason  of  the  inclusion  of  the  steady  forward  speed.  If  the  steady 
forward  motion  is  introduced  to  the  oscillating  ship,  the  possibility  of 
linearization  depends  on  the  hull  shape  parameter  ns  well.  The  main 
difficulty  in  the  oscillating  thin  ship  with  forward  speed  lies  in  the 
fact  that  the  disturbance  generated  by  the  periodical  motion  of  the  ship 
is  weaker  than  the  disturbance  due  to  the  forward  motion.  Consequently, 
the  first  order  theory  would  lead  to  an  unrealistic  conclusion  that  no 

damping  to  the  oscillation  could  exist.  In  order  to  overcome  this 
24 

difficulty,  Newman  employed  two  independent  parameters,  one  of  which  is 
the  oscillation  amplitude  and  the  other  is  the  hull  shape  parameter,  namely 
the  beam-to-length  ratio  of  the  ship.  Although  the  justification  of  the 
damping  and  added  mass  of  the  thin  ship  is  attained  by  the  use  of  two 
parameters,  more  serious  difficulty  appears  when  t he  ship  is  moving  in 
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ambient  waves.  Cross  products  between  the  fluid  velocity  of  incident  waves 

and  that  due  to  the  steady  forward  motion  appear  in  the  same  order  of 

magnitude  as  that  of  the  fluid  motion  due  to  the  oscillation,  making  the 

free  surface  condition  much  more  complicated.  In  this  respect,  the 

25 

application  of  the  slender  body  theory  looks  more  profitable  than  the 

thin  ship  assumption.  The  assumption  of  the  slender  body  increases  the 

order  of  magnitude  of  the  forward  motion  sc  that  it  is  higher  than  that  of 

lateral  or  vertical  motions.  Therefore,  the  effect  of  the  steady  forward 

motion  does  not  appear  at  the  lowest  order  in  the  far  field  expansion  with 

respect  to  the  hull  shape  parameter.  Though  the  effect  of  the  steady 

motion  may  appear  in  the  lowest  order  in  the  expansion  of  the  boundary 

condition  in  the  near  field,  the  first  order  solution  may  take  a  neat  form 

if  the  wave  amplitude  and  the  slenderness  ratio  of  the  ship  are  taken  as 

two  independent  parameters  of  the  perturbation.  However,  some  numerical 

computations  have  revealed  another  difficulty.  This  is  that  the  first 

26 

order  theory  gives  only  an  unrealistic  result.  On  the  other  hand,  it  is 

widely  accepted  that  the  strip  theory  has  been  able  to  present  a  reasonable 
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agreement  with  measurement.  The  strip  theory  is  regarded  in  one  sense 
as  another  slender  body  theory,  although  it  is  originally  derived  by  a 
somewhat  intuitive  method.  In  a  rigorous  sense,  it  is  a  rational  approxi¬ 
mation  for  a  slender  ship  in  oscillations  with  high  frequency  without 
forward  speed.  Thus,  it  is  eventually  known  that  results  may  become 

different  if  the  different  choice  of  magnitude  of  frequency  and  forward 
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velocity  is  taken.  Discussions  in  this  connection  will  be  given  in  the 
next  section. 

Once  the  slender  body  theory  is  employed,  the  boundary  conditions  are 
expanded  by  the  slenderness  ratio.  Terms  of  the  lowest  order  are  taken 
first.  One  of  the  features  of  the  slender  body  theory  is  the  singular 
perturbation.  It  makes  a  difference  in  the  expansion  at  near  field  and 
at  far  field,  and  a  matching  procedure  is  applied  between  them.  The 
problem  which  we  are  going  to  discuss  is  a  slender  body  floating  on  regular 
waves  and  moving  with  a  uniform  average  speed  U  in  the  mean  direction  of 
its  longitudinal  axis.  In  the  most  general  case,  the  direction  of  the 
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forward  velocity  is  different  from  that  of  wave  propagation.  Then  the  ship 
undertakes  oscillations  of  six  degrees  of  freedom  around  its  mean  position, 
among  which  surge,  heave,  and  pitch  are  longitudinal  oscillations  and 
sway,  yaw,  and  roll  are  lateral  oscillations.  Since  the  viscosity  effect 
in  the  longitudinal  oscillations  is  very  little,  the  potential  flow 
theory  seems  to  offer  a  fairly  accurate  prediction,  while  in  the  lateral 
oscillations,  effect  of  viscosity  plays  an  important  role,  so  that  we 
cannot  expect  any  reliable  results  by  theories  without  taking  account  of 
the  viscosity.  Therefore,  we  will  confine  our  discussions  in  the  longi¬ 
tudinal  oscillations  hereafter.  The  general  formulation  for  the  os¬ 
cillation  with  six  degrees  of  freedom  is  described  in  Appendix  A. 

On  developing  the  perturbation  analysis,  we  take  the  ratio  6  of  wave 
amplitude  to  wavelength  and  the  ratio  £  of  beam-to-length  of  the  ship  as 
basic  parameters.  Since  they  are  mutually  independent,  we  can  expand  the 
velocity  potential  with  respect  to  6  first.  The  first  term  is  independent 
of  6  and  represents  the  fluid  motion  when  the  ship  moves  with  uniform 
velocity  in  still  water.  The  linearized  theory  takes  terms  up  to  the  first 
order  with  respect  to  6.  If  we  are  concerned  with  the  ship  motion  in 
regular  waves,  the  term  which  is  linear  to  6  is  a  simple  harmonic  and 
represents  the  oscillatory  part  of  the  velocity  potential.  The  next  stage 
is  the  expansion  of  the  above  portions  of  the  velocity  potential,  which 
have  been  linearized  already  by  6,  by  the  slenderness  ratio  K.  There  is  a 
term  which  is  independent  of  £  in  the  oscillatory  potential.  It  represents 
the  incident  waves  which  may  be  assumed  as  simple  harmonic  too.  The  other 
part  represents  the  disturbance  by  the  ship.  Consider  a  relative  motion 
with  respect  to  the  coordinates  moving  with  the  average  forward  velocity  U, 
and  take  the  axis  of  x  in  the  direction  opposite  to  the  forward  velocity  of 
the  ship  and  the  axis  of  z  vertically  upwards.  Then  the  velocity  potential 
can  be  written  in  the  form  Ux  +  <j>,  and  (j>  satisfies  the  Laplace  equation 

V2  <j>  =  0  (89) 

in  the  space  occupied  by  the  fluid.  The  boundary  conditions  satisfied  by 
the  velocity  potential  are  those  on  the  hull  surface  and  on  the  free 


38 


surface.  If  the  depth  of  water  is  assumed  infinite,  the  condition  of  fluid 
motion  at  infinity,  both  horizontal  and  vertical,  is  that  the  fluid 
velocity  due  to  the  disturbance  by  the  ship  vanishes,  and  that  the  fluid 
motion  is  just  the  sum  of  the  uniform  flow  and  the  incident  wave.  The 
radiation  condition  at  great  distance  should  be  considered  as  well.  If  the 
hull  surface  at  each  instant  is  expressed  by  the  equation 

y  =  f(x,z,t)  (90) 

the  boundary  condition  on  the  hull  surface  becomes 


11  +  (u+  1£\  M  +  M  11  _  i£  =  0 

at  v  3xj  3x  3z  3z  3y  U 


Now  we  consider  only  the  heaving  and  pitching  oscillations  because 

surging  motion  causes  fluid  motion  of  higher  order.  Designate  the  vertical 

displacement  of  the  center  of  gravity  by  z  and  the  angle  of  pitch  by  ^  and 

6 

take  coordinates  x^,  y^,  z^  fixed  to  the  ship.  Then  the  relation  between 
the  coordinate  systems  (x,y,z)  and  (x0,y0,zQ)  is 


x_  =  x  cos  ^  -  (z~z  )  sin  ip 
0  g 

y0  =  y 

z„  =  x  sin  ip  +  (z-z  )  cos  ip 

0  g 


The  equation  of  the  hull  surface  in  reference  to  the  moving  coordinates 
(x0,y0,z0)  is  independent  of  time,  such  as 


yo =  VvV 


It  is  assumed  that  z  and  ip  are  of  the  order  of  6,  and  after  omitting 

g 

higher  order  terms,  one  can  transform  the  hull  boundary  condition  into  the 
following  equation. 
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Now  let  us  examine  the  order  of  magnitude  when  the  ship  is  regarded  very 

slender.  If  n  is  the  slope  of  the  hull  surface  to  the  longitudinal 
x 

axis,  then  its  order  of  magnitude  is  the  slenderness  ratio  C.  We  must  keep 
in  mind  the  fact  that  the  disturbance  velocity  potential  of  a  slender  body 
is  singular  along  its  longitudinal  axis,  and  the  differentiation  of  it  in 
the  direction  of  the  normal  changes  the  order  of  magnitude  by  c  ^  This 
fact  can  be  shown  also  by  adopting  so-called  strained  coordinates  which 
measure  lengthwise  direction  and  lateral  direction  by  different  scales. 

This  procedure  is  well  known  and  will  not  be  repeated  here.  As  a  conse¬ 
quence  of  this  argument,  the  relation  between  the  order  of  magnitude  of 
and  that  of  3n  is 

34>0/9n  =  t_1  0(<J>0) 


40 


2 

Since  n  =  0(e),  the  above  relation  results  <}>  =  0(e  ).  The  oscillatory 

X  U 

part  of  the  fluid  motion  is  composed  of  the  incident  wave  potential  <j)  and 
the  oscillatory  disturbance  potential  The  former  is  independent  of  c 

and  can  be  regarded  as  a  given  function.  In  the  case  of  regular  waves 
propagating  in  the  direction  making  an  angle  a  with  the  x  axis,  the 
velocity  potential  of  the  incident  wave  is  expressed  by 

4>  =  h/g/K  exp[Kz-iK(x  cos  a+y  sin  a)+iwt]  (97) 

w 

where  h  =  wave  amplitude 

K  =  wave  number  2tt/ A 
0)  =  circular  frequency  of  encounter 
The  absolute  frequency  of  the  wave  is 

0)q  =  co  -  U  K  cos  a  (98) 

There  is  a  relation  between  the  wave  number  and  the  frequency  as  follows: 

K  =  Wg/g  =  (w-UK  cos  a)2/g  (99) 

The  order  of  magnitude  of  U  and  w  or  K.  may  not  be  unity,  so  that  we  need  to 

include  these  quantities  in  the  argument  of  the  order.  It  can  be  assumed 

that  the  amplitude  of  the  ship's  oscillation  is  of  the  same  order  as  that 

of  the  wave  amplitude.  The  frequency  of  the  oscillation  is,  however,  not 

the  frequency  of  the  wave  but  the  frequency  of  encounter  to.  The  wave 

amplitude  is  of  the  order  of  6/K  so  that  the  velocity  of  the  oscillatory 

motion  of  the  ship  has  the  order  to  6  K  The  fluid  velocity  of  the 

-1/2 

incident  wave  has,  on  the  other  hand,  the  order  6  K  ,  so  that  they  are 
not  necessarily  the  same  order.  We  have  to  keep  these  facts  in  mind  in 
examining  the  order  of  magnitude  of  the  oscillatory  part  of  the  velocity 
potential.  Taking  the  first  order  terras  with  respect  to  6,  we  obtain 


3fn  9fn  .  3fn  3fn 

4zi-T^(*  -xip)  +  U-r-%  +  U— 

8xq  3zg  '  g  3z0  3x  3z0 


3fn  3  3fo  a  d(t>n 

+  dr  k  <  W  +  jt0  k  (W  ~ u  3z  ^ 


2  2 
3d 


-  JV  <W  +  U<VxW 


\,!0  )z! 


3y3z 


(100) 


If  we  assume  U  and  to  are  both  of  the  order  of  unity,  <J>  will  be  of  the  order 
of  6  £.  The  term  of  the  lowest  order  in  the  above  equation  has  the  order 
of  6  and  the  equation  of  the  same  order  is 


3fn  3fn  3fn  a 

( z  -xjj)  +  U  ^  ip  +  f-  (<p.+<p  ) 
3zq  g  3z„  3z„  3z  T1  ^w 


'0 


'0 


/  3f  32*  32$  \ 

4-  (<p+<p  )  +  u(z  -X*)  I  ~  — -  -^T  J  =  0 

3y  rl  rw'  g  \3z0  9z2  9y3z/ 


(101) 


If  we  write  the  outward  normal  to  the  sectional  form  of  the  hull  in  the 
plane  perpendicular  to  the  x  axis  by  n',  the  above  equation  takes  the  form 


Sn7" 


(z^-xi^-lty) 


3z 
3n  ’ 


U(z^-xt^) 


(102) 


The  right-hand  side  is  regarded  as  a  known  function  provided  that  the 
steady  potential  <t>g  is  known.  One  can  divide  the  periodical  potential 
into  a  part  determined  by  the  ship's  oscillation  and  a  part  originated  by 
the  diffraction  of  the  ambient  wave.  The  former  is  the  radiation  po¬ 
tential  for  which  the  boundary  condition  becomes 
R 
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(103) 


tst  ’  Hr  -  »<«g-*c>  3^  (S2) 

while  the  latter  is  the  diffraction  potential  <p  which  has  to  satisfy  the 
boundary  condition 


^w 

3n'  ~  3n' 


(104) 


Let  us  consider  next  the  boundary  condition  at  the  free  surface.  If 
the  form  of  the  free  surface  is  given  by  the  equation 


z  =  £(x,y,t) 


(105) 


the  kinematical  condition  is 


|i+  (u+|i)|it|*|S.|i  . 

3t  \  3x  /  dx  3y  3y  3z 


(106) 


while  the  condition  of  constant  pressure  is 


Eliminating  t,  between  the  above  equations,  we  obtain 


(107) 


> +  («  f ) 


3_  3£  3_  _3$  1_ 

3x  3y  3y  3z  3z 


If  1(f)  Hi)  I +g: 


008) 


This  is  the  exact  nonlinear  form  of  the  free  surface  condition.  Let  us 
consider  first  the  near  field  and  examine  the  order  of  magnitude  of  each 
term.  We  have  written  the  velocity  potential  in  the  form 


I 


r 


<f>  =  u  4>0  +  4>l  +  *W 


(109) 


We  should  remember  that  the  order  changes  by  e  when  differentiating  <j>  or 
^  by  y  or  z.  Picking  up  the  time- independent  part,  we  obtain 


jV  ,  3 

»  2  9x9y  9y  2  „  2  \  9y  / 


2  3^ 

3*  3y2  u2  32 


+  0(e)  =  0 


(110) 


Then  the  term  of  the  lowest  order  is  (g/U  )  9<(>q/9z  which  has  the  order  of 
e.  Therefore,  the  free  surface  condition  for  the  steady  potential  of  the 
lowest  order  is 


=0  at  z  =  0 


(111) 


that  is  the  condition  of  double  body  flow.  We  cannot  proceed  to  the  next 
step  without  handling  the  nonlinear  terms  in  the  free  surface  condition. 
The  periodical  part  of  the  free  surface  condition  can  be  written  as 


/  9 

9  > 

, 2  3h  , 

„  3f0 , 

/  3  3  \ 

(at  +u 

9x  > 

1  *i  +  “3r  +  2 

"  sr1 

( 9t +u  ur) 

+  0 (6e  )  «  0 


(112) 


If  U  and  to  are  of  the  order  of  unity,  the  lowest  order  term  is  0(d)  and  the 

next  term  is  0(6e).  Although  these  are  linear  with  respect  to  <j>^,  if 

and  <j)  are  assumed  to  be  known  functions,  the  cross  products  between 
w 

derivatives  of  and  4>^  make  the  boundary  value  problem  intractable.  If 
only  the  lowest  order  is  taken,  the  free  surface  condition  for  the  radiation 
potential  becomes 

~  =0  at  z  =  0  (113) 

Thus,  the  double  body  condition  holds  again.  The  boundary  condition  for  the 
diffraction  potential  becomes,  on  the  other  hand, 


8z 


=  -  U 


S2<J V 


*w  .  2 
9z 


(114) 


where  r  is  the  wave  profile  of  the  incident  wave,  since  we  have  the 
relation 


3<}> 

at 


+  u 


9<() 

v 

9x 


+  8  5. 


(115) 


Now  let  us  construct  the  boundary  value  problem.  The  field  equation 
is  the  Laplace  equation,  but  the  governing  equation  in  the  near  field  is 
reduced  to  the  two-dimensional  form  in  the  y-z  plane  because  of  the  singular 
perturbation.  Therefore,  the  boundary  value  problem  is  to  find  a  plane 
harmonic  function  with  the  boundary  condition  at  the  hull  surface,  where  the 
normal  velocity  is  prescribed,  and  with  the  condition  9<j)^/9z  =  0  at  z  =  0 
for  the  radiation  potential.  The  condition  of  infinity  i;  left  unspecified, 
yielding  an  indefiniteness  to  the  solution.  Now  we  assume  the  solution  of 
the  two-dimensional  prob1 em  of  the  Laplace  equation 

2  2 

t±  +  = 

2  2 
9y  dzZ 


0 


(116) 


r . 


with  the  boundary  condition  on  the  hull  surface.  We  have  omitted  the  term 

2 

0(£  )  in  reducing  the  Laplace  equation  to  the  two-dimensional  form.  Then 
we  can  include  the  term  0(e)  which  may  be  the  second  order  in  the  pertur¬ 
bation  expansion.  Next  the  inner  solution  for  the  slender  body  can  be 
expressed  by 


4>n  =  <p^2D)  +  g1(x)  +  z  g2(x)  +  y  g3(x)  (117) 

where  Reference  20,  is  the  two-dimensional  solution  of  the  boundary 

value  problem,  and  g^(x),  g2(x),  and  g^(x)  represent  the  indefiniteness  of 

the  solution  which  may  be  determined  by  the  matching  procedure  with  the  far 

field  behavior  of  the  fluid  motion.  In  the  above  solution,  we  include 

2 

terms  up  to  the  second  order,  but  we  have  omitted  terms  of  0 (e  )  in  the 
expansion  of  the  free  surface  condition,  so  that  the  third  and  fourth  terms 
on  the  right-hand  side  of  the  above  expression  must  be  deleted.  Then  we 
are  obliged  to  take  the  expression 


<t>N  =  4>^2D^  +  gx(x)  (118) 

The  free  surface  condition  in  the  far  field,  on  the  other  ha. id,  takes  a 
different  form  because  the  differentiation  with  respect  to  y  or  z  does  not 
affect  the  order  of  magnitude.  Therefore,  the  leading  term  of  the  equation 
gives  simply  the  linearized  boundary  condition  as  follows. 

(I?+Ul;)2  «■  +  eff  ■  0  atz-°  <119> 


For  the  periodical  motion  with  circular  frequency  w,  we  can  write  the 
complex  form  as 


I 


and  the  boundary  condition  on  the  free  surface  becomes 


3$  3*. 

V2iU“3T+U  sr=  0 

dx 


The  Green’s  function  for  this  boundary  condition  is 


G(x,y,z;x' ,y' ,z')  =  7-  +  ^ 


rl  r2 


00  c 

ij  “» f 


dn  exp[-jy-y’|/  m^+n^  -  i  m(x-x')] 


_oo  0 


•  {cos(nz+e)  cos(nz'+e)  -  cos  nz  cos  nz'}//  in  +n 


^  j  dm  exp[(z+z')  (mU-Ko)  ^/g- 1  y-y '  |  J m^-(mU-Ho)^/g^ 


(121) 


-im(x-x’)]  (mU +w>)^  /  / m^-(mU+u))^/g^ 


(122) 


where 


^  =  /  (x-x,)2+(y-y')2+(z-z,)‘ 


r  =  /  (x-x*)  +(y-y ' )  +(z+z') 


(123) 


tan  e  =  -  (raU-co)  /gn 


If  the  radical  in  the  last  integral  becomes  imaginary,  it  has  an  appropriate 
sign  in  accordance  with  the  radiation  condition.  The  outer  solution  is 
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expressed  by  means  of  this  Green's  function.  This  can  be  achieved  by  the 
application  of  Green's  theorem  in  similar  fashion  to  the  problem  of  steady 
forward  motion  discussed  in  the  former  section.  Because  of  the  slender 
body  assumption,  the  first  order  solution  takes  the  form  in  which  the 
singularities  on  the  hull  surface  shrink  to  a  line  distribution  of  wave 
sources  along  the  x-axis.  Then  the  far  field  potential  is  expressed  by 

4>p  =  -  e1U)t  Jm(x')  G(x,y,z;x'  ,0,0)  dx'  (124) 

where  m(x)  is  the  source  density  which  is  determined  by  the  matching 
procedure  between  the  near  field  and  far  field  solutions.  Thanks  to  the 
simple  condition  30/9z  =  0  at  z  =  0,  the  source  density  can  be  determined 
in  a  simple  way  such  as 

m(x)  =  ~  \t\  (iuH"u  (zg-x^-Cw)]  (125) 


where  B(x)  is  the  width  at  the  water  plane  at  each  transverse  section  of 
the  ship.  The  unknown  function  g-^(x)  is  determined  by  the  inner  expansion 
of  the  far  field  potential  and  the  final  result  becomes 


-  *<2D>  -  * 


sgn(x-x')  £n(2|x-x'|)  dx' 


f  d  m(x' 

J  dx' 

J*m(x')  G' (x, 0, 0;xf , 0, 0)  dx ' 


(12b) 


where 


G' (x,y,z;x' ,y' ,z')  =  G(x, y,z;x ' , y ' ,z ' )  -  - -  —  (127) 

rl  r  2 

We  have  assumed  that  m(x)  vanishes  at  both  ends  of  the  ship.  One  can 
formulate  the  forces  and  moments  acting  on  the  ship  by  the  integration  of 
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pressure  on  the  hull  surface,  such  as  the  added  mass  and  damping.  Numerical 
computations  are  not  so  simple  and  may  require  much  computer  time,  if  the 
forward  velocity  is  taken  into  account,  but  they  are  possible  at  any  rate. 
There  has  been  some  attempt  of  calculating  the  added  mass  and  damping  of 
a  ship  making  forced  heaving  and  pitching  oscillations  during  forward 
motion  in  still  water.  However,  it  has  failed  to  obtain  any  useful  result 
because  the  numerical  results  show  very  unrealistic  features,  in  spite  of 
a  rational  appearance  of  the  formulation.  This  fact  will  be  discussed  in 
another  section. 

Various  Cases  of  the  Order  of  Magnitude 
of  Frequency  Parameter  and  Froude 
Number 

In  the  preceding  section,  we  have  assumed  that  the  frequency  parameter 
of  oscillations  and  Froude  number  are  both  of  the  order  of  unity,  and  it 
was  expected  that  the  theory  was  valid  without  any  restriction  in  magni¬ 
tudes  of  frequency  and  forward  speed.  However,  the  results  were  quite 
disappointing.  It  was  noted,  on  the  other  hand,  that  the  different  choice 
in  the  assumption  of  the  order  of  magnitude  of  frequency  parameter  or 
Froude  number  might  result  in  different  formulation.  The  discussion  of  the 
order  of  magnitude  in  the  boundary  condition  may  derive  different  solutions 
for  different  assumptions.  Here  we  will  consider  cases  that  the  Froude 
number  is  not  so  large,  or  the  frequency  of  oscillations  is  not  so  small. 

Now  let  us  begin  with  the  case  of  low  Froude  number.  The  Froude 

number  of  conventional  merchantile  vessels  is  not  much  greater  than  0.3, 

2  -1  -1 

so  that  the  speed  parameter  U  /g£  =  Yq  is  of  the  order  of  10  ,  which  may 

be  regarded  as  a  small  parameter.  Therefore,  let  us  assume  that  U/i^£  = 

1/2 

0(c  ).  The  order  of  magnitude  of  the  radiation  potential  is  6c  but  the 

1/2 

effect  of  the  forward  speed  appears  in  the  term  of  the  order  6  c  .  The 
lowest  order  term  in  the  free  surface  condition  lias  the  order  of  6  and  the 
next  order  is  6c.  If  we  take  up  to  the  order  of  6c,  the  free  surface 
condition  in  the  near  field  becomes 

324>  34>  3  2<f> 

7y  +  «8r  +  Uf-w7T=o  at2  =  0  (128> 

dt  d  Z 


I 
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Therefore,  we  can  take  account  of  the  effect  of  forward  speed  in  the 
boundary  condition  on  the  hull  surface  without  bothering  about  the  quadratic 
terms  in  the  free  surface  condition  making  the  boundary  value  problem  in¬ 
tractable.  The  near  field  potential  for  the  radiation  potential  is 
determined  from  the  two-dimensional  boundary  value  problem  with  the 
linearized  free  surface  condition  for  an  oscillating  body  such  as 


3z 


(129) 


The  well  known  method  of  solution  for  the  two-dimensional  problem  can  be 
applied  to  the  determination  of  the  two-dimensional  potential  The 

solution,  however,  is  not  identical  to  that  for  an  oscillating  cylinder  of 
infinite  length  because  of  the  term  involving  the  forward  speed.  The 
results  from  the  boundary  condition  on  the  hull  surface, 

3^  -  Uu<Vx*)-U*)  -  »<V*«  (jpj  (130) 

The  last  term  on  the  right-hand  side  may  add  some  complication.  The 
velocity  potential  in  the  near  field  is  then  expressed  by 

<t>N  =  4>^20^  +  g^x)  +  z  g2(x)  (131) 


The  term  yg^(x)  which  appeared  in  Equation  (117)  is  omitted  because  of  the 
symmetry  of  longitudinal  oscillations.  It  is  readily  shown  that  there  is 
a  relation 


S2(X)  =  gl(x) 


(132) 


because  of  the  free  surface  condition.  The  function  g^(x)  is  determined  by 
matching  with  the  inner  expansion  of  the  far  field  potential.  Since  we 
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have  assumed  that  the  Froude  number  is  of  the  order  of  e  ,  we  expand 
the  far  field  potential  which  has  been  given  in  the  preceding  section  by  U 
and  discard  terms  of  higher  order.  The  result  is 


g1(x)  =  -  2 j  dm,^i"  sgn(x-x')  £n(2|x-x'|)  dx' 

+  7Tvj*m(x’)  IHq(vIx-x'  !  )-Yq(v|x-x'  | )  ]  dx' 

+  2  tt  i  nj— H0(v|x-x’|)-Y0(v|x-x'|) 

+v|x-x’|  -Hl(v|x-x'  IHY^vIx-x'  |)j  J  dx' 

+  2  TT  i  v(l+4ft2)J  m(x')  e2ivfi(x~x')  H^2)(v|x-x'  |)  dx' 

-  4  tt  v  ftfm(x’)  e21^(x-x')  HJ2)(v|x-x’  |)  dx'  ( 


(133) 


where  v 


ft  =  gj  U/g 

Hq  =  Struve  function 

(2) 

=  Bessel  functions  of  the  second  and  third  kinds 

Because  of  the  term  zg„(x),  the  boundary  value  problem  becomes  a  little 

2  (2D) 

different  from  the  pure  two-dimensional  solution  for  <J)  .  The  most 

remarkable  feature  of  this  case  is  that  the  two-dimensional  solution  is 

related  to  the  free  surface  condition 


V  *(2W  -  D 


(134) 


instead  of  that  for  the  double  body  flow  given  in  the  preceding  section. 
This  is  the  consequence  of  taking  the  second  term  of  the  perturbation 
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expansion.  Though  the  free  surface  condition  does  not  involve  the  forward 

velocity,  the  hull  boundary  condition  does. 

Another  choice  is  the  high  frequency  case.  The  order  of  magnitude  of 

-1/2 

the  frequency  parameter  is  assumed  to  be  C  .  In  the  first  place,  we 
assume  the  Froude  number  still  remains  in  the  order  of  unity.  Because  of 
the  relation 


V  =  2tt  _  £ _ 

K  \  2  2 

2  U  cos  a 


(l+2fi  cos  (x-  /l+4\2  cos  a) 


(135) 


the  ratio  of  wavelength  to  the  ship's  length  is  of  the  order  of  t  ",  so 
that  the  wave  is  not  extremely  short.  Since  the  wave  slope  is  of  the 

1/9 

order  of  6,  the  wave  amplitude  is  of  the  order  of  6  t  ".  The  velocity  of 
the  orbital  motion  of  wave  is  0(6),  so  that  the  order  of  is  6  o.  In  the 
boundary  condition  on  the  hull  surface  for  the  radiation  potential,  namely 


a<t>R  .  ,  ..  3z  ...  3z  ...  ..3  /3*0 

3iT  =  1  3^  -  u  *  3^  -  U(  V  3T 


the  first  term  is  of  the  order  of  6,  while  other  terms  which  are  related  to 

1/2 

the  forward  speed  are  of  the  order  of  6  c  .  On  the  other  hand,  the 
lowest  order  of  the  free  surface  condition  is  6  and  the  next  is  6  c^".  If 
we  take  the  lowest  order  terms  only  in  both  boundary  conditions,  the  so¬ 
lution  is  in  the  case  of  zero  forward  speed.  We  can  discuss  the  effect  of 

1/2 

the  forward  speed  by  taking  terms  of  order  6c  .  In  this  case,  the  free 
surface  condition  for  the  radiation  potential  becomes 


3^  3(f)  ?2<pn  3ch 

+  2  U  3t3x  +  2  U  37'  3t3y  "  U  ~ 7  3t  =  0 


Because  of  the  terms  except  the  first  and  second  ones,  we  cannot  solve  the 

boundary  value  problem  in  two  dimensions  by  ordinary  methods.  Ogilvie  and 
29 

Tuck  employed  a  successive  method  by  which  the  solution  for  the  above 


boundary  condition  was  derived  from  the  solution  for  zero  forward  speed. 

The  solution  obtained  involves  integrals  over  the  plane  z  =  0.  It  seems 
of  some  interest  that  the  outer  solution  for  the  radiation  potential  takes 
a  degenerated  form  like 

<f>R  ~  |m(x)+2ift(z-ijy[)  dm^  j  exp  [  iiot+v(z-i  j  y  |  )  ]  (138) 

The  first  term  in  parentheses  means  a  simple  harmonic  plane  wave  propagat¬ 
ing  outwards  in  the  y  direction,  while  the  second  term  is  related  to  the 
variation  of  the  source  density  along  the  x-axis  which  means  the  inter¬ 
ference  between  different  sections.  An  extensive  discussion  of  this  case 
is  given  by  Ogilvie  and  Tuck  and  will  not  be  reproduced  here. 

Next  we  consider  the  case  of  high  frequency  with  low  forward  speed, 

-1/2  1/2 

that  is  to  =  0(e  ),  U  =  0(e  ).  This  is  the  case  of  short  waves  such 

that  the  wavelength  is  of  the  same  order  as  that  of  the  breadth  of  the  ship. 

3/2 

In  this  case,  the  order  of  magnitude  of  ^  is  6  e  .  The  free  surface 
condition  becomes 


a2<f>,  a<j>, 

+  g  ^  =  0  (139) 

dtZ  az 

1/2  3/2 

up  to  the  order  of  <5  e  .  The  next  term  is  of  the  order  of  6  £  ,  which 

includes  the  effect  of  the  forward  speed,  but  has  quadratic  forms.  There¬ 
fore,  the  ordinary  linearized  treatment  can  be  applied  to  the  first  order 
only,  by  which  the  effect  of  the  forward  speed  cannot  be  taken  into 
account.  The  asymptotic  form  for  the  outer  solution  for  the  radiation 
potential  takes  the  form 


— eiWt  +  V(z~ilyl)  m(x) 
g 


( 140) 


that  means  outward-going  plane  waves.  There  are  no  three-dimensional  terms 
involved  and  the  solution  is  purely  two-dimensional.  The  strip  theory  Is 
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y 


then  valid  but  without  any  effect  of  the  forward  speed.  It  should  be  noted 
that  the  diffraction  problem  cannot  be  treated  in  the  same  way,  because  the 
assumption  of  the  slow  variation  along  the  x-axis  is  no  longer  valid.  The 
diffraction  in  short  waves  will  be  discussed  in  another  section. 


Hydrodynamic  Forces  in  Heaving  and 
Pitching 

It  is  widely  known  that  the  hydrodynamic  forces  and  moments  acting  on 
oscillating  ships  are  predicted  by  strip  theory  with  fairly  good  accuracy. 
However,  the  reliability  of  the  strip  theory  is  still  open  to  doubt  as 
discrepancies  between  computed  and  measured  results  are  observed  occa¬ 
sionally.  These  discrepancies  may  be  attributed  to  the  effect  of  the 
forward  speed  and  fluid  motion  in  three  dimensions.  In  the  preceding 
section,  we  have  observed  that  both  the  forward  speed  and  the  three- 

dimensionality  can  be  taken  into  account  within  a  plausible  approximation 

1  /2 

of  the  perturbation  scheme,  if  the  forward  speed  is  of  the  order  of  c 
where  e  is  the  beam-to-length  ratio. 

The  hydrodynamic  forces  are  obtained  by  the  integration  of  pressure 
over  the  hull  surface.  The  fluid  pressure  in  the  near  field  is  given  by 


-  p  (p-Po}  = 


i  0) 


+  U 


34>] 

9x 


+  U 


30 

A 

3y 


(2D) 


3$ 

w 


(2D) 


+  L! 


30 

dT 


(2D) 


30 


(2D) 


(141) 


up  to  the  order  of  6c.  The  hydrostatic  pressure  is  omitted  because  it  is 
simply  determined  by  the  geometrical  relations.  Although  the  third  and 
fourth  terms  are  omitted  in  the  usual  linearized  theory  for  oscillating 
ships,  they  appear  in  the  same  order  of  magnitude  as  the  second  term,  so 
that  they  have  to  be  retained  if  one  wishes  to  take  account  of  the  effect 
of  the  forward  speed.  The  vertical  force  is  given  by  the  integral 


F 


2 


(P-P0}  In  dS 


(142) 
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Although  the  quadratic  terms  in  Equation  (141)  seem  to  be  troublesome  in 

30 

evaluating  the  above  integral,  a  theorem  which  has  been  proven  by  Tuck 
becomes  a  powerful  aid.  The  theorem  is  proved  in  Appendix  B.  Consider  a 
velocity  vector  V  of  an  irrotational  motion  of  an  inviscid  fluid  outside 
the  hull,  which  satisfies  the  boundary  condition  on  the  hull  surface  such 


V  •  n  =  0 


(143) 


where  n  is  a  unit  vector  along  the  outward  normal  to  the  hull  surface. 
Next  we  define  the  vector  m  by  the  relation 


(144) 


Jc  * 

Further  we  define  vectors  n  and  m  by 


n  =  r  x  n 


2*  =  -  h  (£xn) 


(145) 


where  r  is  the  position  vector  (x,y,z).  Then,  the  following  relations  are 


valid. 


[nxHrKV*  V<f>)  ]  dS 


=  -  n  <P  w 


11 


*<p+n* (V'Vtp)  ]  dS 


f  * 

=  -  I  n  <j>  w 


(146) 


where  w  is  the  z-component  of  V  and  Lp  is  the  still  waterline  of  the  hull. 
The  line  integral  on  the  right-hand  side  is  omitted  when  the  body  is 
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slender  and  wall-sided  at  the  water  plane.  For  the  vertical  force,  we  take 
the  z-component  of  the  first  equation,  while  for  the  pitching  moment,  the 
y-component  of  the  second  equation  is  taken.  In  the  present  case,  we  put 


V 


<  *0  3(*>o 

,u’  u  IT  *  u  3^ 


(147) 


*  * 

Next  we  define  vector  functions  x>  tP»  X  >  and  V  ,  which  are  two-dimensional 
harmonic  functions  in  the  lower  half  space  outside  the  ship  and  satisfy  the 
boundary  conditions  on  the  hull  surface,  such  as 


3X 

dip 

9n  =  -’ 

U  -r—  =  m 
dn 

* 

* 

3X  * 

M  * 

3^=  2  * 

U  3^  =  2 

(148) 


Furthermore,  these  functions  are  assumed  to  satisfy  the  free  surface 
condition 


a  - v  *  =  ° 


at  z 


(149) 


We  have  expressed  the  near-field  expression  of  the  radiation  potential  for 
a  slender  ship  by 


<t>N  =  4)('2D^  +  (l+vz)  gx(x) 


(150) 


Then,  the  vertical  force  is  written  in  the  form 

.(2D)  3z 


Fz  "  P 


(  dx 

*  J  ( 


1  0) 


9n 


r  dc 


C(x) 


f  dx 

f 

u 

. 

ajt.<2D)  ^  ^2D>  a4<2Di  ^  »*o2W  n*20’ 


dx 


dy  dy 


dz 


dz 


dc  (151) 
t)n  (coni.) 


+  dx  {iu)g1(x)+Ug^(x) }  j 
J  '  C(x) 


(1+vz)  |^r  dc 


(151) 


where  C(x)  is  the  contour  of  each  section  along  which  the  integral  with 
respect  to  c  is  taken.  On  applying  the  above  mentioned  theorem  to  the 
second  term,  we  obtain 


F  = 
z 


,  ,  (2D)  3z  , 

1  <*>  <t>  g^T  dc 


=  p  f  dx  f 
J  "u(x) 

+  pfdxf  U3^(a^)* 

J  JC(x) 


+  p  J{itog1(x)+Ug|(x)}  {-B(x)+vS(x)  }  dx 


(2D) 


dc 


(152) 


The  last  term  is  derived  by  the  application  of  Gauss'  theorem  in  which 

B(x)  and  S(x)  are  breadth  at  the  waterline  and  sectional  area  of  each 

transveri*  section,  respectively.  If  we  write  the  z-components  of  x  and  £ 

by  x  and  ijj  ,  respectively,  the  two-dimensional  part  of  the  radiation 
z  z 

potential  is  expressed  as 


<j/2D^  =  {iu>(z  -xifO-Uip}  x,  +  U(z  -x\p)  p 

6  Z  g  Z 


(153) 


Then,  the  vertical  force  is  written  as 


F  =  -  p  I  dx  {w2(z  -xij>)+iwlty}  Xz  dc 

J  8  JC(x) 

r  r  9*z 

+  i  p  co  U  I  dx  (z  ~xp)  I  p^  dc 

^  8  u(x) 


dp 


r  r 

-  i  p  0)  U  dx  (zg-x\Ji)  Xz  g^r 

J  JC(x) 

-  p  U2Jdx  (zg-4)  f  ^  ^  dc 

JC(x) 
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dc 


(154) 
(cont . ) 


% 


(154) 


{ia)g1(x)+Ug^(x)  }  {-B (x)+vS (x)  }  dx 


Because  of  Green's  reciprocal  relation,  which  can  apply  to  x^  and  \p  ,  such 
as 


3^ 

dc 


the  second  and  third  terms  are  cancelled.  Therefore,  we  obtain 


“Pl 


dx  {co  (z  -xi|j)+ia)UiJj} 
g 


U2J'dx  (zg-xt|0  J" 


9^z, 

Xz  3n'  dc 


C(x) 


-  p  U  |  dx  (zg-xt p)  |  \pz  dc 

C(x) 


+  P 


-iw{B(x)~VS(x) }  +U  ^  (B(x)-vS(x). 


e,  (x)  dx 


1/2 

The  last  term  is  derived  by  integration  by  parts.  If  we  assume  U  =  0(t  ) 

2 

and  omit  the  term  of  0(6e  ),  the  second  term  drops  out.  A  similar  ex¬ 
pression  is  obtained  for  the  hydrodynamic  moment  about  the  y  axis.  The 
first  term  indicates  the  result  obtained  by  the  strip  theory,  and  the  other 
terms  give  the  effect  of  the  forward  speed  and  the  three-dimensional  motion 


Numerical  Results  of  Radiation  Problem 
at  Zero  Forward  Speed 

Although  numerical  analysis  works  when  the  forward  speed  and  the 

effect  of  three-dimensionality  are  present,  no  published  result  of  this 

kind  is  known  so  far.  There  is  a  rather  comprehensive  result,  on  the  other 

hand,  for  the  case  of  zero  forward  speed  by  means  of  a  similar  formu- 
31 

lation.  It  is  obtained  simply  by  letting  U  =  0  in  the  original  formula 
for  finite  speed  with  no  substantial  difference  in  the  method  of  numerical 
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calculation.  The  result  can  well  illustrate  the  effect  of  three- 


dimensionality  of  the  fluid  motion  which  is  remarkable  at  lower  frequencies. 

In  the  first  place,  let  us  put  (p  =  e1Wt  $  and  express  the  two- 
dimensional  solution  for  a  heaving  cylinder  in  the  form  as 


.  (2D) 

$  =  -  a. 
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kz  cos  ky  .  Kz 

e  —  dk-nie  cos  Ky 
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(2m-l)  .' 
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~2m-l  /  \  ~2m-2 

8  z  \+KS  — ' 


„  2m-l  2  2 

dz  \z  +y 


„  2m-2  2  2 

dz  \  z  +y 


(156) 


It  is  readily  shown  that  this  expression  satisfies  the  free  surface 
condition 


-T--  -  K  <J>  =  0  at  z  =  0 

dZ 


(157) 


2 

where  K  =  V  =  oj  /g  in  the  present  case.  Since  the  inner  expansion  of  the 
above  function  is 


*<2D>  = 


[£nKr+y+Kr  cos  0 


(l-£nKr-y)+Kr0  sin  0+tt i ] 


+ 


cos  2m0+ 


Kr 
2m- 1 


cos 


(2m-l)  0 


(158) 


where  we  have  employed  the  cylindrical  coordinates 

z  =  -  r  cos  0,  y  =  r  sin  0  (159) 

and  y  is  Euler's  constant,  0.5772157. 

In  order  to  make  the  inner  expansion  of  the  far  field  potential  match  the 
above,  we  employ  the  expression  in  the  near  field  as 


$  =  $(2D)  _  a0(l+Kz)  (y+TTi) 


-"2  ( 1+Kz )  3q(x' )  sgn(x-x')  £n(2K|x-x'|)  dx' 


Ji 

|  K(l+Kz)  I  aQ(x’)  [H0(K|x-x'|)+Y0(K|x-x'|) 

»  n 


+2iJQ(K|x-x' | ) ]  dx' 


(160) 


where  x  =  +£  is  the  x  coordinate  at  each  end  of  the  ship.  The  second  term 
on  the  right-hand  side  is  added  to  conform  to  the  expansion  of  the  two- 
dimensional  solution.  We  can  rewrite  the  above  expression  as 


$  =  $(2D)  + 


y  (1+Kz)  I  a^(x)  N(K|x-x'|)  sgn(x-x')  dx" 


(161) 


where 


fu  ru 

N(u)  =  -y-£n2u+-“  I  Hq(u')  du'  +  j  YQ(u')  du'  "  71  1 


+  it  i  Jq(u')  du' 


(162) 


It  can  be  shown  that 


£im  N(K | x-x' | )  =  0 
K-k» 


The  coefficient  Aq  can  be  determined  by  the  method  well  known  in  the  two- 
dimensional  theory  of  a  heaving  cylinder.  Then  we  can  determine  a^  by 
solving  the  above  equation.  However,  we  need  not  solve  the  equation 
exactly.  The  reason  is  as  follows.  If  the  frequency  is  very  low  or  to  the 
contrary  very  high,  the  integral  on  the  right-hand  side  can  be  neglected. 

So  we  can  put 


aQ(x)  =  V (x)  Aq(x) 


(168) 


It  may  be  assumed  that  the  deviation  from  the  above  relation  at  inter¬ 
mediate  frequencies  is  not  large.  Therefore,  we  approximate 


a0(x> 


£ 

V(x)-  i  K  I  {V' (x' ) AqCx' )+V(x* )Aq(x* ) }  N(K | x-x' | ) sgn(x-x' )dx' 

-it 


x  Ao(x^ 


(169) 


In  order  to  determine  other  coefficients,  we  substitute  V(x)A(x)  for  aQ  in 
the  boundary  condition  and  put 


£ 

W(x)  =  ~  K  (V (x' ) AQ(x' )+V(x' )Aq(x' ) }  N(K | x-x' | ) sgn(x-x' )dx'  (170) 

-£ 


Then  the  boundary  condition  becomes 


30 

9n ' 


(2D) 


=  [V(x)-W(x) ] 


3z 
3n ' 


(171) 


The  solution  of  the  two-dimensional  problem  with  this  boundary  condition 

determines  other  coefficients  a„  .  We  can  rewrite  the  function  ari(x)  and 

zm  U 

W(x)  in  the  following  form  for  numerical  purposes. 
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aQ(x)  =  V (x)  Aq(x)  1-  y  KAq(x) {N(k| £+x | )+N(k| £-x | ) } 


J6 

-  ~  K2  A0(x)  I  [V(xf ) Aq (x* )-V (x)Aq(x) ]  N'(K|x-x' |)dx' 

»  n 


W(x)  =  |  K  V(x)  AQ(x)  {N(K{  5l+x|  ) -HN (K 1 5,-x  [  )  } 


(172) 


f*2j 


[V(x’ )Aq(x' )-V(x) Aq(x) ]  N'(K[x-x' [ )dx' 


(173) 


where 


N’(K|x-x'|)  =  ~  Tx-xTf  +  fH0(K[x-x'j)  +  j  Y0(K|x-x'|) 

+  TTi  Jq(K[x-X*  |) 

The  vertical  component  of  the  force  acting  on  the  ship  by  the  fluid 
pressure  is  given  by 


(174) 


'-//It 


(175) 


and  the  moment  about  the  y-axis  is 


My  =  JJ  (It  *-  It  Z)  P  dS 


(176) 


The  second  term  in  the  parentheses  can  be  omitted  because  of  the  slender 
body  assumption.  Then  the  vertical  force  and  the  pitching  moment  on  the 
slender  ship  are  given  by 


L 
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ds  = 


x  dx 


(177) 


where  C(x)  is  the  contour  of  each  transverse  section.  The  function  f (x) 
gives  the  force  per  unit  length  at  the  section.  If  we  divide  the  velocity 
potential  in  two-dimensional  and  three-dimensional  portions  as 


♦  -  e1"*  (*(2W+*(3D>) 


(178) 


we  can  write 


f(x)  =  i  p  u)  j  0^2D^  ds  +  i  p 

JC(x) 


to 


.(3D)  3z  , 

*  8^  ds 


C(x) 


=  fL(x)  +  f  2 (x) 


(179) 


Now  take  the  added  mass  m^  and  the  damping  coefficient  for  a  heaving 
cylinder  of  infinite  length.  Then  the  two-dimensional  portion  of  the 
sectional  force  is  expressed  by 


f. (x)  =  -  (ion  +N  )  [V(x)-W(x)] 
1  z  z 


(180) 


while  the  three-dimensional  portion  has  the  expression 


f„(x)  =  i  p  to  S (x)  -  i  £&  B(x)  W(x) 

z  to 


(181) 


where  S(x)  is  the  area  and  B(x)  is  the  waterline  width  of  each  section. 
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In  pure  heaving,  we  put  V(x)  =  iu)z^,  and  we  can  write 


f  (x)e 


itot 


-  (m  +m')  ( — to  z  )  -  (N  +N')  i  to  z 
z  z  g  z  z  g 


(182) 


where  m'  and  N'  indicate  the  three-dimensional  effect  to  the  added  mass 
z  z 

and  damping  coefficient,  respectively.  In  the  case  of  pure  pitching,  we 
put  V(x)  =  -  itdxiji  and  we  can  write 


f  (x)e 


iwt 


-  (m  x+m")  to  +  (N  x+N")  i  w 
z  z  z  z 


(183) 


Here  we  put 


a  =  m  +  m  b  =  N  +  N 

Z  Z  Z  Z 


*  * 
d  =  m  x+m"  e  =  N  x  +  N" 
z  z  z  z 


(184) 


and  define  the  following  integrals. 
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(183) 
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These  constants  give  hydrodynamic  coefficients  in  the  coupled  equation  of 
heaving  and  pitching  as  follows. 


(a+pV)  V  +  b  z  +  c  z  -  dV  -  e\p  -  g\p  =  F 

b  b  O 

(A+k2  pV)  V  +  BiJ>  +  OJ>  -  D  *z*  -  E  z  -  G  z  =  M 
yy  8  g  8 


(186) 


where  V  is  the  displacement  volume  and  k  is  the  radius  of  gyration.  On 

yy 

account  of  Haskind's  relation,  we  have  the  relation 


D  =  d  E  =  e 


(187) 


For  numerical  example,  the  hull  form  of  Series  60,  C  =  0.7  is  employed 

D 

because  reliable  data  of  model  experiments  have  been  available  for 
comparison. 

*  *  *  * 

The  coefficients  a  ,  b  ,  d  ,  e  ,  and  also  a,  b.  A,  B,  d,  and  e  are 
calculated.  They  are  compared  with  results  computed  by  means  of  the  strip 
theory  and  the  measured  results,  as  shown  in  Figures  5  through  7.  One  can 


Figure  5  -  Hydrodynamic  Coefficients  a  and  b 


observe  a  remarkable  improvement  in  agreement  with  measured  results  by  the 
calculation  including  the  three-dimensional  effect,  especially  at  lower 
frequencies.  One  of  the  reasons  of  deviation  with  the  results  by  the 
strip  theory  may  be  attributed  to  the  logarithmic  term,  S-nKr,  in  the  two- 
dimensional  solution  which  makes  the  added  mass  infinite  at  zero  frequency 
limit.  The  three-dimensional  part  of  the  velocity  potential  has  a  term 
which  cancels  the  above  logarithmic  singularity. 

The  application  of  the  simple  slender  body  theory,  which  has  been 
described  in  an  earlier  section,  presents  such  an  unrealistic  result  as 
negative  value  of  added  mass  and  infinite  increase  in  damping  coefficient 
at  higher  frequencies.  Computed  results  for  the  source  term  a^Cx)  are 
shown  in  Figure  8.  The  result  by  the  simple  slender  body  theory  which 
determines  the  near-field  solution  by  the  condition  for  a  double  body 
shows  much  deviation  from  other  theories.  This  fact  may  be  the  main  reason 
of  the  ill  behavior  of  the  slender  body  calculation.  The  fair  agreement 
between  strip  theory  and  the  present  calculation  suggests  that  the  source 
term  can  be  determined  by  the  two-dimensional  calculation  without  regard 
for  the  three-dimensional  effect,  or  aQ(x)  =  V(x)Aq(x). 


Figure  8  -  Source  Distribution 


As  explained  in  the  preceding  section,  the  present  theory  is  a  conse¬ 
quence  of  the  inclusion  of  the  second  term  in  the  slender  body  expansion. 
Nevertheless  it  is  by  no  means  a  second  order  theory,  but  still  a  first 
order  theory.  The  formal  perturbation  procedure,  taking  successive 
approximations,  starting  from  the  lowest  order  term  and  proceeding  to  the 
second  order  approximation,  never  leads  to  the  same  result  as  the  above. 
This  fact  is  a  peculiar  feature  of  the  asymptotic  expansion  of  the  singular 
solution  such  as  the  present  problem. 

WAVE  PRESSURE  ON  SLENDER  SHIPS 

Boundary  Value  Problem  for  the 
Diffraction  Potential 

Although  the  strip  theory  is  employed  in  the  usual  practice  of  pre¬ 
dicting  wave  exciting  forces  on  ships,  the  diffraction  problem  does  not 
admit  the  use  of  the  strip  theory  in  longitudinal  waves.  Although  the 
strip  theory  is  an  acceptable  approximation  in  high  frequencies  for  the 
radiation  problem,  the  short  wavelength  associated  with  the  high  frequency 
invalidates  the  condition  of  slow  variation  along  the  ship's  axis.  In  the 
case  of  long  waves,  on  the  other  hand,  the  frequency  becomes  low  and  the 
three-dimensional  effect  comes  in  as  in  the  radiation  problem. 

Now  we  consider  first  the  case  of  a  ship  with  forward  speed  in  long 
waves,  so  that  we  assume  U  and  w  are  both  of  order  of  unity.  The  boundary 
value  problem  for  the  diffraction  potential  has  been  given  previously,  but 
its  solution  needs  some  contrivance.  The  boundary  Value  problem  in  the 
near  field  is  the  two-dimensional  Laplace  equation 
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with  the  hull  surface  boundary  condition  of  Equation  (104) 


!!? 

an '  9n ' 


(188) 


(189) 
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and  the  free  surface  condition  of  Equation  (114) 


% 

~  =  -  U  ^  ^  at  z  =  0  (190) 

dz  w  ^  2 

dz 


If  the  incident  wave  propagates  along  the  direction  of  x,  £  is  constant  in 

w 

the  Lransverse  plane.  Then  the  function 


*'  =  *0 


+  U( 


'w  3z 


(191) 


is  a  plane  harmonic  function  in  the  lower  half  space.  The  boundary 
condition  for  it  is 


3<j> 

Tv 

3n’ 


+  U 


(192) 


on  the  hull  surface,  and 


3z 


0 


at  z  =  0 


(393) 


If  we  take  only  the  first  order  terms,  we  can  write  the  hull  surface 
boundary  condition  as 


3$' 
3n ' 
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3z 

3n' 


+  U  L 


3 

3n ' 


(194) 


or 
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(195) 
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This  is  the  same  form  as  the  hull  boundary  condition  of  the  radiation 
potential  for  heaving,  because  the  latter  is 
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3z 
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(196) 


Therefore,  the  diffraction  is  taken  into  account  in  the  hull  boundary 

condition  by  replacing  the  vertical  movement  of  the  section  by  the  relative 

displacement  to  the  surface  of  the  incident  wave.  This  relation  holds  in 

1/2 

the  case  of  U  =  0(e  )  too.  Then  we  can  take  up  to  the  next  term,  so  that 

the  relation  is  to  be  modified  as 


a<j>' 

8n' 


-  i  ai  £  (1+Kz)  ~~r  +  u  C  A 
w  8n  w  9n' 


(197) 


Kz 

The  added  term  Kz  comes  from  the  exponential  factor  e  in  the  incident 
wave  potential  and  indicates  the  Smith  correction. 

The  boundary  value  problem  for  the  diffraction  potential  is  now  re¬ 
duced  to  a  similar  form  to  that  for  the  radiation  problem.  The  field 
equation  for  the  potential  <(>'  is  the  two-dimensional  Laplace  equation 


i±'  +  fr.  =  o 

3y2  az2 


(198) 


with  the  hull  boundary  condition  given  above.  The  free  surface  condition 
at  z  =  0  is 


a<j>' 

dz 


0 


when  U  =  0(1),  u)  =  0(1) 


or 
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when  U  = 


0(e1/2), 


(199) 


oj  =  0(1) 


t 

Therefore,  the  inner  solution  for  the  diffraction  potential  in  the  case  of 
u  =  0(1)  is  expressed  by  the  form 


*o  ■  »'l2W  - "  +  (1+r z) 800  (200) 

where  cp'  tjle  golution  0f  the  two-dimensional  problem  with  the  above 
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mentioned  boundary  conditions.  The  factor  to  z/g  in  the  third  term  on  the 

1/2 

right-hand  side  can  be  acftied  only  when  U  =  0(e  )  and  should  be  omitted 
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in  the  case  of  U  =  0(1).  Another  case  to  be  considered  is  co  =  0(e  ). 

This  is  the  short  wave  case,  but  some  complication  appears  if  we  apply  the 

slender  body  theory.  The  basic  idea  of  the  slender  body  is  that  the  field 

equation  in  the  near  field  can  be  reduced  to  the  two-dimensional  Laplace 

equation.  However,  the  short  wavelength  hinders  the  above  possibility.  In 
-1/2 

the  case  of  cj  =  0(e  )  with  U  =  0(1),  the  ratio  of  the  wavelength  to  the 

1/2 

ship's  length  is  =  0(e  )  and  the  variation  of  the  flow  field  along 

the  x-axis  is  related  to  the  wavelength.  If  the  order  of  the  diffraction 
potential  is  6  £,  the  order  of  magnitude  of  each  term  in  the  Laplace 
equation  in  the  near  field  is 


32^d 

3  < 

3x2 

3y2 

; 

3z‘ 

(6) 

(6e_1) 

(<5e 

Therefore,  omitting  the  term  of  higher  order  than  £,  we  get  a  two- 

dimensional  Laplace  equation.  However,  the  omitted  term  in  the  long  wave 

2 

case  has  been  of  higher  order,  £  .  Therefore,  the  validity  of  the  two- 

dimensional  equation  becomes  much  weaker  in  comparison  with  the  long  wave 
* 

case.  This  may  damage  the  accuracy  appreciably.  The  free  surface 
condition  in  this  case  is 
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(201) 
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if  terms  of  order  6  e  are  retained.  The  terms  involving  the  effect  of 
steady  forward  potential  prevents  the  straightforward  solution,  so  that 

this  case  should  not  be  used  for  the  practical  purpose  of  prediction. 

1/2 

If  the  forward  speed  is  low,  namely  U  =  0(e  ),  the  ratio  of  wave 

length  to  ship's  length  is  of  the  order  of  e.  Since  each  term  in  the 
three-dimensional  Laplace  equation  has  the  same  order  of  magnitude,  its 
two-dimensional  version  is  no  longer  valid.  Therefore,  the  strip  theory  is 
not  applicable  to  the  diffraction  problem  in  the  longitudinal  waves.  An 
alternative  method  for  the  diffraction  problem  in  short  waves  will  be 
discussed  later. 


Wave  Pressure  and  Hydrodynamic  Forces 

As  was  mentioned  before,  the  diffraction  problem  requires  not  only  the 
integrated  total  force,  but  some  local  quantities  such  as  wave  pressure  at 
each  point  on  the  hull  surface  and  the  distribution  of  forces  along  the  x- 
axis.  If  we  write 


=  <J>  _  +  4> 

D  Tw 


(202) 


the  periodical  pressure  on  the  hull  surf^e  is  gi^en  by 


I(p.Po).lu*  +  „|i  +  „|o|i  +  ul?0|i 


(203) 


up  to  the  order  of  <5  e,  if  U  and  to  are  both  of  the  order  of  unity.  Al¬ 
though  the  third  and  fourth  terms  are  omitted  in  the  usual  theory,  they 
appear  in  the  same  order  as  the  first  and  second  terms,  so  they  have  to  be 
taken  into  account  if  the  effect  of  forward  speed  is  considered.  The 
vertical  force  is  given  in  the  same  way  as  in  the  radiation  problem,  so 
Tuck's  theorem  can  be  applied  again.  The  vertical  force  is  given  by 
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At  a  great  distance  from  the  ship,  the  diffraction  potential  and  the 
incident  wave  potential  <j>  satisfy  the  linearized  free  surface  condition 


+  g 


84) 

9z 


at  z 


(205) 


However,  the  boundary  condition  for  at  z  =  0  changes  near  the  ship  as 
was  pointed  out  before.  Instead  of  4>D>  we  can  define  a  function  <p'  which 
satisfies  the  above  free  surface  condition  even  in  the  near  field.  We 
assume  that  4>p  coincides  with  4>D  in  the  far  field,  while  it  coincides  with 
<p’  in  the  near  field.  Then  in  the  near  field,  we  can  put 


*  =  *D  +  *w 

=  K  ~ u  ar  +  *w  (206) 


Since  the  boundary  condition  on  the  hull  surface  is 


3£ 

3n 


0 


the  boundary  condition  for  <p^  becomes 

Hr  „  r  <l(Hr)  H* 

3n  4w  3n  \3z  /  3n 


(207) 


Next  we  assume  the  auxiliary  functions  Xz  and  satisfy  the  boundary 
condition 


at  z  =  0 


(208) 


and  evaluate  the  integrals 
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We  have  the  relation  which  is  derived  by  Green's  identity. 
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where  Z  is  a  vertical  cylinder  of  large  radius  surrounding  the  ship.  We 
have  omitted  a  line  integral  term  which  can  be  assumed  small  because  of 
the  slender  body  assumption.  Because  of  the  boundary  condition  on  the  hull 
surface,  there  is  a  relation 


JJ  xz  k  ds  -  JJ  u Is  (S2 


Xz  ds 


while  the  radiation  condition  satisfied  by  <p^  and  Xz  gives 
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Therefore,  we  have  the  relation 
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Therefore,  also 
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The  same  relations  hold  when  x  is  replaced  by  \p  .  Then  the  vertical  force 

z  z 

is  expressed  by 
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If  there  is  no  forward  speed,  we  have 


F  =  i  p  u) 
z 


(211) 


This  relation  is  called  the  Haskind  relation.  The  functions  y  and  i b  are 

z  z 

related  to  the  solution  of  the  radiation  problem.  Therefore,  the  above 
relation  enables  the  wave  excitation  to  be  evaluated  without  solving  the 
diffraction  problem.  Similar  relations  are  obtained  for  pitching  moment 
too. 
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Wave  Pressure  in  Short  Waves 

-  —1/2 

When  the  order  of  magnitude  of  the  frequency  parameter  w/Jl/g  is  e  , 

the  variation  of  the  diffraction  potential  along  the  x-axis  is  no  longer 
small.  Therefore,  the  slender  body  assumption,  in  which  the  flow  field 
varies  slowly  along  the  axis  of  the  body,  does  not  hold,  and  the  strip 
theory  is  by  no  means  applicable.  We  need  not  worry  about  this  fact,  as 
far  as  the  total  force  or  moment  is  concerned,  because  the  total  force  be¬ 
comes  very  small  and  unimportant  in  short  waves.  However,  we  need  another 
theory  if  we  intend  to  discuss  more  local  phenomena  such  as  pressure  on  the 
surface  of  the  body. 

In  the  case  of  a  long  ship  in  a  heading  with  the  propagation  of  short 

waves,  one  may  think  of  an  infinitely  long  cylinder  placed  in  waves  with 

its  axis  parallel  to  the  propagation  of  the  wave.  However,  it  has  been 
32 

shown  by  Ursell  that  there  is  no  steady  state  solution  in  this  case. 

The  fluid  motion  is  highly  three-dimensional  and  cannot  be  replaced  by  the 

two-dimensional  solution.  However,  there  is  an  acceptable  assumption  which 

can  reduce  the  boundary  value  problem  to  a  much  more  simplified  form. 

That  is,  the  variation  of  the  diffraction  potential  in  the  longitudinal 

axis  of  the  ship  deviates  from  the  sinusoidal  variation  very  slightly. 

-1/2 

Here,  we  consider  the  case  of  oj  =  0(e  ). 

As  mentioned  before,  it  is  not  easy  to  formulate  the  linearized  so- 

-1/2. 

lution  of  the  boundary  value  problem  in  the  case  of  w  =  0(e  )  and 

1/2 

U  =  0(1).  If  U  =  0(e  ),  the  effect  of  the  forward  speed  appears  only  in 

the  higher  order  term.  Therefore,  we  will  confine  our  discussion  in  the 
case  without  forward  speed.  The  velocity  potential  of  the  incident  wave 
which  propagates  in  the  positive  x-direction  can  be  expressed  by 


=  c  h  exp(Kz-iKx+iiot) 


(212) 


where  c  =  g/w  is  the  phase  velocity  of  the  wave.  Now  let  us  express  the 
diffraction  potential  in  the  form 


<j>n  =  ip(x,y,z)  exp(i'jjt-iKx) 


(213) 


According  to  the  assumption  stated  above,  the  function  ip(x,y,z)  varies 
very  slowly  in  the  x-direction.  In  order  to  facilitate  the  analysis,  we 
take  the  Fourier  transform  of  4*^  such  as 
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where  \p  is  the  Fourier  transform  of  i p.  The  diffraction  potential  satisfies 
the  Laplace  equation 
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and  its  Fourier  transform  is 
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Substituting  the  above  mentioned  expression  for  the  diffraction  potential, 
Equation  (214),  shifting  the  first  term  of  Equation  (216)  to  the  right- 
hand  side  and  putting  k  -  K  =  k',  we  obtain 


ff+  (k.+K)2? 

3y  3z^ 


(217) 


The  basic  assumption  of  the  slow  variation  of  suggests  that,  if  \p  «*  0(1) 
then  <Kk',y,z)  =  0(l/k').  Therefore,  k'  must  be  0(1)  if  =  0(1).  Since 
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il 
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K  =  0(e  )  and  the  singular  perturbation  in  the  near  field  suggests 

3/ 3y  =  0(e  ^) ,  3/3z  =  0(e  ^),  we  obtain  for  the  lowest  order,  the 
Helmholtz  equation 


^f+^f=K2? 

3y  3z 

or,  by  the  Fourier  inversion, 
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(218) 


(219) 


The  omitted  terms  are  of  the  order  higher  by  e.  If  we  employ  the  cylindri¬ 
cal  coordinates  r,  0, 


z  =  -  r  cos  0,  y  =  r  sin  0 


(220) 


particular  solutions  of  the  differential  equation  are 

I  (Kr)  \  i  cos  n0 

n  (221) 

K  (Kr)  1  l  sin  n0 
n 

where  1^  and  are  Bessel  functions  of  imaginary  arguments. 

Since  the  forward  velocity  is  absent  in  the  boundary  value  problem,  the 
boundary  condition  on  the  free  surface  for  the  diffraction  potential 
applies  to  the  function  ip  too. 


dip 

3z 
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at  z  =  0 


(222) 


One  can  construct  a  solution  which  satisfies  the  above  equation  and 
vanishes  a*:  infinity  by  combining  the  basic  solutions  of  Equation  (221) 
as  follows. 
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(223) 


V  =  K.  „(Kr)  cos(2n-2)  8  +  2  K.  .  (Kr)  cos(2n-l)  6 
n  zn-z  zn-l 

+  K„  (Kr)  cos  2  n  8 
2n 

However,  this  type  of  function  does  not  form  a  complete  set  and  we  cannot 

express  the  diffraction  potential  which  satisfies  the  body  boundary 

condition  by  means  of  the  linear  combination  of  these  functions.  In  order 

to  find  the  missing  term,  we  consider  a  distribution  of  wave  sources  with 
icot mm  i  Kx 

density  1/2  o(x)  e  along  the  x-axis.  Taking  the  Fourier  transform 

of  o(x)  such  as 
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we  have  the  velocity  potential  of  the  type  e  where 


. ,,  Kz  f  exp  [  iax 
X  =  -  xKe  — US— 

-K  / 


.  /2  2 

-l/K  -a 


ill  a, 


/F-c? 


a(-a-K)  da 


+  K  e 


Kz 


J  -J1 

K 


exp [ lax-/a2-K2 1 y ] ]  ~(_x_ 

/  2  2 
/a  -K 


K)  da 


+  ^  |  eiax  ‘o(-a-K)da 


C  ™  /  2  2  J~2  2  .  J  2  ~  2 

I  /u  -a  cos/u  -a  z+K  sin/ii  -a 

J  .2  “2  "2  "  1 


-u  y 


du 


u  -  a  +  K 


(225) 


Now  we  change  the  variables  of  integration  to 


,  „  /  2  2 
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and  assume  the  slow  variation  of  c(x),  such  as 


o(k)  =  0(l/k) 

and  Ky  =  0(1),  Kz  =  0(1),  because  the  wavelength  is  comparable  with  the 
breadth  of  the  body.  Then  the  asymptotic  expression  of  the  above  integrals 
for  large  K  leads  to  the  following  expression. 
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where  x  =  -  l  is  the  forward  end  of  the  body.  The  second  term  is  the  two- 

dimensional  part  of  the  potential,  while  the  first  term  gives  the  three- 

dimensional  effect.  Because  of  the  presence  of  the  factor  ^  in  the  first 

1/2 

term,  the  order  of  the  second  term  is  higher  by  C  .If  we  take  the 
lowest  order  term  only,  we  can  omit  the  second  term.  Then  the  boundary 
value  problem  becomes  extremely  simple.  The  boundary  condition  on  the 
hull  surface  is 
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the  boundary  condition  can  be  written  as 
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It  is  known  that  the  above  condition  can  be  satisfied  by  putting 


i  iKx  . 

^  =  e  x  (231) 

and  the  density  of  the  wave  source  is  determined  by  such  a  simple  equation 
as 


/2nK  e"l7rM  |  ^ 

J-h  ^ 


This  is  Abel's  integral  equation,  the  solution  of  which  is 

,  iTT/4 

o(x)  =  -  (233) 

Att3K(x+Z) 
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This  approximation  was  actually  given  by  Faltinsen  who  calculated  the 

pressure  and  force  acting  on  a  body  with  semicircular  cross  section. 

However,  the  solution  shows  the  density  of  the  source  being  infinite  at 

the  forward  end  of  the  body  which  results  in  an  infinite  pressure  there. 

Such  things  never  happen  in  actual  phenomena.  In  order  to  eliminate  this 

difficulty,  one  has  to  retain  other  terms  in  the  expression  of  ip  other 

than  the  lowest  order.  In  this  case,  the  function  ip  is  expressed  by  a 

linear  combination  of  x  and  T  as 

A  n 

v  =  Ak  eKz_i"/4  f  Z&-  dr 

_oo 

(oo 

V1  \  eKZ  +  JL  pn  \)  (234) 

n=  ] 
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where 


-  2  7T  K I  y  I  eKz  (235) 

The  coefficients  can  be  determined  independently  from  o(x)  in  the 

boundary  condition,  while  a(x)  should  satisfy  an  integral  equation  of  the 

Volterra  type  of  the  second  kind.  The  solution  of  this  integral  equation 

never  presents  a  singular  behavior.  Numerical  results  for  the  pressure 

distribution  on  a  body  with  semicircular  cross  section  are  compared  with 

measurements  as  well  as  the  results  of  the  lowest  order  approximation  and 
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those  by  the  strip  theory.  Fairly  good  agreement  with  measured  results 
is  obtained  by  the  present  calculation  as  shown  in  Figures  9  and  10,  while 


Figure  9  -  Distribution  of  Wave  Pressure 
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Figure  10  -  Wave  Pressure  at  Square  Station  Number  8 

remarkable  discrepancies  are  observed  in  the  results  by  the  lowest  order 
approximation  and  those  by  the  strip  theory.  It  should  be  noted  that  the 
present  solution  can  never  be  derived  by  the  successive  method  starting 
from  the  lowest  order  approximation.  This  then  is  another  example  which 
illustrates  the  importance  of  the  inclusion  of  terms  next  to  the  lowest 
order  in  the  asymptotic  expansion. 

ADDED  RESISTANCE  IN  WAVES 

Formulation  of  the  Steady  Drift  Force 
in  Waves 

When  a  ship  floating  on  water  encounters  incident  waves,  the  hydro- 
dynamic  pressure  on  the  hull  surface  is  periodical  which  yields  periodical 
forces  and  moments  acting  horizontally  on  the  ship.  They  are  called  wave 
excitations  or  wave  loads. 
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However,  the  ship  experiences  another  kind  of  horizontal  force  which 
is  independent  of  time.  Such  a  force  is  usually  called  the  steady  drift 
force  by  the  wave.  When  the  ship  is  travelling  among  ambient  waves,  the 
steady  drift  force  results  in  an  increase  of  resistance  to  the  forward 
motion  of  the  ship.  This  is  called  the  added  resistance  in  waves.  In  some 
cases,  the  added  resistance  attains  an  amount  even  larger  than  the  re¬ 
sistance  when  the  ship  travels  in  calm  water.  It  is  generally  regarded 
that  the  added  resistance  should  occupy  a  major  part  of  the  sea  margin 
which  should  be  considered  when  the  power  of  ships  in  the  service  condition 
is  predicted  from  model  test  results.  This  problem  looks  to  be  more 
complicated  when  compared  with  the  problem  of  resistance  in  calm  water  or 
periodical  hydrodynamic  forces  in  oscillations,  because  the  steady  drift 
force  originates  from  the  periodical  action  of  waves.  The  model  test  in 
waves  is  difficult  also  in  respect  to  achieving  accurate  measurements.  It 
needs  much  time  and  labor.  A  reliable  method  of  prediction  by  theoretical 
formula  has  been  expected  for  a  long  time.  From  the  theoretical  side,  this 
problem  looks  difficult  because  the  steady  drift  force  belongs  to  the 
second  order  forces.  However,  the  problem  of  added  resistance  is  now 
regarded  as  one  of  the  most  successful  applications  of  hydrodynamic  theories 
in  the  practical  field  of  ship  building. 

Since  the  drift  force  is  a  kind  of  second  order  phenomena,  one  may 
think  of  the  necessity  of  the  second  order  solution  of  the  boundary  value 
problem,  but  it  is  not  the  case.  The  steady  force  can  be  calculated  from 
the  pressure  integral  on  the  hull  surface.  In  doing  so,  however,  we  have 
to  pick  up  terms  of  the  second  order  without  exception.  This  is  quite  a 
difficult  task,  because  we  need  solutions  up  to  the  second  order  complete¬ 
ly.  Fortunately  we  have  another  method  which  does  not  need  the  second 
order  solutions.  That  is  the  momentum  analysis  by  assuming  a  large  surface 
surrounding  the  ship  at  a  great  distance  as  a  momentum  control  surface. 

The  steady  horizontal  force  is  evaluated  by  the  time  average  of  the 
momentum  flux  across  the  control  surface.  Its  principle  can  be  easily 
understood  by  considering  a  two-dimensional  case.  Here,  we  consider  a  cy¬ 
lindrical  body  of  uniform  cross  section  which  is  floating  on  the  surface  of 


water  and  a  train  of  regular  waves  comes  in  the  direction  perpendicular  to 
the  axis  of  the  body.  The  incident  wave  is  partly  transmitted  beyond  the 
body  but  partly  reflected,  generating  reflected  waves  which  propagate  in 
the  direction  opposite  to  the  propagation  of  the  incident  wave.  At  a 
great  distance  from  the  body,  there  exist  regular  sinusoidal  waves.  On  the 
weather  side,  there  are  incident  waves  and  reflected  waves,  while  on  the 
lee  side,  there  are  transmitted  waves.  If  we  assume  vertical  planes  on 
both  sides  at  a  great  distance  from  the  body,  the  momentum  flux  across 
these  planes  can  be  evaluated  simply  by  the  expression  for  regular  waves. 

If  we  write  the  amplitude  of  the  incident  wave  by  h,  that  of  reflected 
wave  by  h^,  and  that  of  transmitted  wave  by  h^,,  there  is  a  relation  for  the 
energy  conservation  law 


h2  =  h^  +  h2  (236) 

Now  we  apply  the  momentum  principle  to  the  fluid  between  the  vertical 

surfaces  and  take  the  time  average  for  one  period  of  the  wave.  The  result 

shows  the  steady  horizontal  force  D  in  the  direction  of  the  propagation  of 

the  incident  wave  experienced  by  the  floating  body,  the  amount  of  which  is 
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given  by  a  simple  relation 


D  =  f  P  g  h2  (237) 

Since  the  first  order  forces  are  periodic,  the  steady  force  is  the  second 
order  force.  Nevertheless,  it  is  calculated  by  the  first  order  solution 
of  the  diffraction  problem.  A  similar  analysis  can  be  applied  to  the 
three-dimensional  problem  even  when  the  forward  speed  is  present. 

Instead  of  the  actual  case  in  which  the  ship  penetrates  into  waves 
with  a  velocity  U,  we  consider  a  ship  floating  on  a  uniform  stream  with  a 
train  of  regular  waves.  A  constant  horizontal  force  is  assumed  in  such  a 
way  that  keeps  the  average  position  of  the  center  of  gravity  of  the  ship 
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at  a  fixed  point.  The  velocity  potential  of  the  fluid  motion  in  this  case 
is  expressed  by  Ux  +  <p  as  before.  Consider  a  vertical  circular  cylinder 
of  large  radius  with  its  axis  through  the  origin  of  the  coordinates.  Tak¬ 
ing  the  immersed  part  of  the  surface  of  the  ship  S,  the  cylindrical 
surface  E,  and  the  portion  of  the  free  surface  E^,  inside  the  surface  E,  as 
the  momentum  control  surfaces,  we  apply  the  momentum  principle  to  the 
fluid  contained  in  the  space  bounded  by  these  surfaces  and  consider  the 
fact  that  there  is  no  flux  across  S  and  E^  and  the  pressure  is  constant 
on  Eq.  Then  the  rate  of  change  of  momentum  M  of  the  fluid  enclosed  by 
these  surfaces  is 


dM 

dt 


n  dS  + 


n  dS  +  p 


V  $(n*V4>)  dS 


+  (gravitational  force) 


(238) 


where  n  is  the  unit  normal  to  each  surface  directing  inward  to  the  fluid 
under  consideration,  and  we  have  put 

$  =  Ux  +  (}>  (239) 


The  force  acting  on  the  ship  is  the  integral  pressure  on  S.  Taking  the  x- 
component,  the  force  in  the  x-direction  is  given  by 


F 

x 


p  n  dS 
r  x 


(u+  (U  nx+n*V<J>)  dS  (240) 

E  "E 

where  subscript  x  means  the  x-component  of  each  vector. 

If  the  time  average  is  taken,  dfWdt  vanishes  because  of  the  periodicity 


dt 


p  n  dS  +  p 
r  x 
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(U  nx+n*V$)  dS  =  0 


where  the  bar  means  the  time-averaged  during  one  period. 

The  average  force  in  the  x-direction  yields  the  resistance  to  the  forward 
motion  of  the  ship. 


D  =  F 


•  JJ  -  % 


dS  +  P  fa  (u  nx+2*^<f>)  dS 

•'•'e 


Taking  cylindrical  coordinates  R  and  0, 


x  =  R  cos  0,  y  =  R  sin  0 


making  use  of  the  pressure  equation 


p  p0  1  . ..  2  lT 
-  =  _  _x  _  _  (V<f>)  -  U  ^  -  gz 


3t  2 


and  designating  the  elevation  of  the  free  surface  by  we  obtain 


D=M  Rd6J  +  i 

J  0  -00 


y  (V<£)  +gZ  j  COS  0  dz 


Because  of  the  periodic  motion,  we  can  put 


2tt  0 

f  R  d0  J  +gz )  cos  0  dz  =  0 


Jhen  we  have  the  relation 


2ir 

-pj 

J0 

R  d0 

J 

C  [iili  _ 

|  1  3x  3R 

-00 

I  (V<J>)2  cos  0  dz 

2tt 

If  COS  6 

2t[  _ 

+  P  f 

R  d0 

dz  -  pg  j  C2  R  cos  0  d0 

(246) 

•'o 

J 

0 

"  J  0 

where  the  free  surface  elevation  becomes 


5 


ifM 

g  Lat 


+  \  <V4>)2 


+u 


3£ 

3x 


(247) 


Thus,  the  lowest  order  terms  are  quadratic  functions  of  the  disturbance 
velocity  potential.  A  similar  formulation  is  applied  to  the  y-component 
of  the  steady  force,  giving 


F 

y 


2tt  £ 

*"  {  [fl-j 


|9i  1 

1  at  2 


(V<())  +gz  sin  0 


dz 


(248) 


On  developing  the  formulation,  a  relation  of  the  energy  conservation  is 
utilized  as  well.  The  energy  contained  in  the  space  bounded  by  the  control 
surface  is 


E  = 


+gz 


dV 


(249) 


If  the  control  surface  is  moving  with  normal  velocity  vn  di  ected  inward, 
the  rate  of  change  of  the  energy  is  given  by 


clE  =  £  fff  JL 

dt  2  J  J  J  at 


(V4>r  dV 


+gz 


v  dS 
n 


(250) 


Since  we  have 


ijjj  it  <w>2  dv-  /{]«*>  (,f)dv 


by  Green's  theorem  and 


1  /T7A  %2  ^  3$  1  ,  . 

2  (V$)  +  gz  =  9F  "  p  ^-Pq5 


from  the  pressure  equation,  one  can  write 


dE  f  f  f  3<J>  3$  ,  /  3$  p  \ 

3t  J  J^3n  3t  \  3t  p  )  vn 


dS 


30 

3n 


on  S  and  by  the  boundary  condition,  while  vn  =  0  on  E  and  p  =  0  on  EQ. 
Then  the  rate  of  change  of  energy  becomes 


dE 

dt 


p  v  dS 
r  n 


30  30  ,e 

3ti  ¥t  dS 


(251) 


If  the  time  average  is  taken,  the  first  term  on  the  right-hand  side  gives 
the  work  done  by  the  ship.  The  ship  is  floating  freely  on  the  free  sur¬ 
face  and  no  external  force  exists  except  the  constant  towing  force  and  the 
gravitational  force  which  keep  the  average  position  of  the  ship  fixed  in 
space.  Therefore,  no  work  is  done  nor  is  there  any  dissipation  of  energy, 
because  the  viscosity  is  neglected.  Therefore, 


dS 


0 


(252) 
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Owing  to  the  periodicity  of  the  motion 


f 

i 

\ 

I 


dE 

dt 


0 


(253) 


we  have  the  relation 


|»|£dS  =  0 
3n  3t 


or 


(254) 


The  velocity  potential  is  the  sum  of  the  incident  wave,  the  steady  forward 
motion,  and  the  periodical  disturbances,  so  we  can  write 

<}i  =  4>w  +  U  +  <J>i  (255) 


Now  we  consider  that  the  direction  of  the  wave  propagation  which  makes  an 
angle  a  with  the  direction  of  the  forward  motion  of  the  ship,  namely  the 
negative  direction  of  x.  The  incident  wave  potential  is 

<J>  =  c  h  exp(Kz+iKx  cos  a+iKy  sin  a+iu>t)  (256) 

w 


The  contribution  of  <J>q  to  D  is  the  wave  resistance  in  calm  water.  The 
periodical  disturbance  potential  has,  on  the  other  hand,  an  asymptotic 
expression 


2  i  e 


iu)t 


-  8-7T/2 

0+7T/2- 

-  -7T/2 

tt/2  - 

a1  exp[a..z+ia..R  cos(a-0)] 
H(a^ ,a)  — - - - - -  da 


/u4?r 


cos  a 


(257) 
(cont . ) 
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e+*n 

2  i  eiu)t  1  H(a2,a) 

Vir/2 


a2  exp[a2z+ia2R  cos(a-0)] 


/l^fT 


da 


cos  a 


+  0(R-1) 


(257) 


where 


al 

a2 


g  1  -  2  (2  cos  a  ±  A  -4(2  cos  a 

2  2 
U  2  cos  a 


(258) 


The  interval  -  <  a  <  (Xq  must  be  omitted  from  the  integral  where 

*  0  for  _<_  1/ 4 
=  cos  (1/4(2)  for  (2  >  1/4 


The  function  H(a±,a)  is  the  Kochin  function  which  is  determined  by  singu¬ 
larities  representing  the  hull.  Inserting  the  above  expressions  for  the 
velocity  potential  in  the  integral  giving  the  force  in  x-direction,  and 
taking  the  time  average  of  the  lowest  order  term,  we  obtain  an  expression 
for  the  mean  resistance  to  the  forward  motion.  It  is  composed  of  the  wave 
resistance  in  the  calm  water  and  the  added  resistance  due  to  incident 
waves.  Let  us  define  the  latter  by  AR.  The  derivation  of  the  final 
formula  is  rather  lengthy  and  will  not  be  reproduced  here.  Only  the  result 
is  shown  as  follows. ^ 


AR  =  2  u  p 


3.(8,  cos  0-K  cos  a) 

— — --  ■  -  -  d0 

/l-4(2  cos  0 


+  2 


.p  J 

a. 


2ir-a, 


|H(a,,0) 


_  a0(a0  cos  0-K  cos  a) 

i2  2  2,  - — 

/l-4(2  cos  0 


(259) 
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where  a  is  replaced  by  6  in  and  a2*  Note  that  this  formula  is  valid 
when  the  energy  relation  mentioned  before  holds.  If  there  is  another  kind 
of  energy  dissipation  such  as  viscosity,  the  above  relation  needs  some 
alteration.  In  order  to  facilitate  numerical  work,  the  integration 
variable  is  changed  to 


m  =  K, 


1  -  2  cos  6  ±  /1-4&  cos  8 
2  cos  6 


(260) 


where  Kn  =  g/U.  Then  we  obtain 


AR  =  4 


r  -K-  , 

r  r  (m+K.ar 

-  J  +J  —r= 

L  -°°  K2j  /(m+KQl 


(m-K  cos  a) 


|H*(m)|2  dm  (261) 


v2  2 
Q)  -KQm 


where 


H  (m)  =  H(ai,0)  i  *  1,  2 


K1  “  I  V1+2J»^>’  K2  =  \  K0(l+2^-/l+4n) 


(262) 


(263) 


If  ft  <  1/4,  the  interval 


K0(l-2ft-/f=4ft)  <  m  <  j  K  (l-2ft+/L-4fl) 


must  be  omitted  from  the  integral. 

Simplification  of  the  Formula 

At  a  great  distance  from  the  ship,  the  disturbance  by  the  ship,  in¬ 
cluding  both  radiation  and  diffraction,  is  represented  by  a  combination  of 
wave  sources  and  y-directed  wave  dipoles  distributed  along  the  x-axis. 
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The  former  represents  the  longitudinal  disturbance  such  as  heaving  and 
pitching,  and  the  latter  represents  the  lateral  disturbance  such  as  sway¬ 
ing,  yawing,  and  rolling.  Now  let  us  designate  densities  of  these  wave 
sources  and  wave  dipoles  by  m(x)e  and  Ji(x)e^wt.  Then  the  Kochin 
function  is  expressed  as 


J1 
{m(x)+iai  sin  0y(x)}  e 


i  a.  x  cos  6 


dx 


(264) 


Then  we  put 


f  1  ai 

H^a^e)  =J  m(x)  e 

f  1  ai 

H2(ai,0)  -I  p(x)  e 


i  a.  x  cos  0 


dx 


i  a.  x  cos  0 


dx 


(265) 


On  integrating  with  respect  to  0  the  Integrals  in  Equation  (264),  odd 

o 

functions  of  0  vanish.  Therefore,  we  can  replace  |H(ai>0) |  by 


|H1(ai,0)|2  +  a2  sin2  0  |H2(a±,0) | 


Changing  the  integration  variable  from  0  to  m,  we  obtain 


AR  =  4  ir  p 


-K, 


■J  */ 

-»  -k2  J 


(m+K-fi)  (m-K  cos  a)  *  „ 

- —  — -  |H  (m)r  dm  (266) 

/(m+K_f2)4-K2m2 


where 


|  H*  (m) 


(m+K  (2)4  -  K2  m2  *  , 

jH1(m)|2+ - - - 2 - ^ -  |H2(m)|2 


(267) 


and 
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H^(m)  ® J  m(x)  e*  m  X  dx 
H*(m)  *  J" p(x)  e1  m  x  dx 


(268) 


One  can  evaluate  the  above  Integral  provided  that  the  densities  of  sources 
and  dipoles  are  known. 

Considerable  data  of  computation  of  added  resistance  In  head  sea  waves 
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by  means  of  the  above  formula  have  been  reported.  As  for  the  case  of 
oblique  seas,  numerical  results  so  far  obtained  are  rather  infrequent  be¬ 
cause  of  the  complication  of  ship  motions.  Unfortunately  the  above  formula 
has  failed  to  strike  the  fancy  of  engineers  of  a  practical  mind,  because 
the  integral  in  the  infinite  interval  necessitates  much  calculation.  It 

will  be  found  that  the  formula  can  be  simplified  if  high  frequencies  such 

—1/2  38  1/2 

as  co  =  0(e  )  are  assumed.  If  we  assume  first  U  =  0(e  '  )  and 

-1 II  -1 

co  =  0(e  '  ),  then  =  0(e  )  and  0  =  0(1).  Omitting  terms  of  0(e),  the 

formula  of  Equation  (266)  is  reduced  to 


"p  J 


(m-K  cos  a) 


g 


|H*(m) 


(269) 


Substitute  Equation  (268)  in  the  above  and  apply  the  Fourier  integral 
theorem.  Now  we  assume  that  m(x)  and  y(x)  vanish  at  both  ends  of  the  ship. 
Then,  by  integration  by  parts,  we  have 


Jm(x)  eimxmdx*-ij" 


dm(x)  i  m  x  , 
— e  dx 

dx 


The  added  resistance  due  to  the  longitudinal  disturbance  is  then  repre¬ 
sented  by 


OO 

=  A  IT  p  I  (m-K  cos  a)|H^(m)|^  dm 
.~00 

-  -  8  tt2  p  Re  j*  |^K  cos  a|m(x)  1 2+im(x)  dm ^  j  dx 


(270) 


W 


where  m(x)  is  the  complex  conjugate  of  m(x).  Here  we  write 

arg  m(x)  =  61(x) 


(271) 


Then  we  obtain 


AR 


1 


8  n 


K  cos  a+ 


dfi^x) 

dx 


dx 


(272) 


A  similar  reduction  is  applied  to  the  portion  by  the  lateral  disturbance. 
In  consequence,  the  added  resistance  is  given  by 


AR  =  -  8  7T 


2 


2 


K  cos  a+ 


dfi^(x)  | 
dx  I 


+  |u(x) |2 

g 


K  cos  a+ 


d62(x)  '■ 
dx  j 


dx 


(273) 


where 


arg  p(x)  =  62(x)  (274) 

It  can  be  readily  shown,  that  the  expression  for  AR  is  valid  in  the  case 

-1/2  1 

of  U  =  0(1)  and  u'  =  0(e  )  too,  but  the  expression  for  the  added 

resistance  due  to  lateral  disturbances  takes  another  form  which  is  more 
complicated  and  spoils  the  practical  utility  to  some  extent.  A  similar 
simplification  can  be  applied  to  the  lateral  drift  force,  and  the  final 
result  is 


D 

y 


8  ir  p 


1 


r  4 

2^j  |m(x) | |y(x) |sin{51(x)-62(x) } 
-  g 


-K  sin  a 


dx 
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(275'* 


* 


The  densities  of  sources  and  dipoles  are  determined  from  the  inner  so¬ 
lution  at  the  near  field  which  has  been  discussed  previously.  They  are 
identical  with  the  source  term  or  dipole  term  in  the  solution  for  two- 
dimensional  cylinders  in  the  corresponding  motion  of  each  transverse 
section.  However,  the  existence  of  the  forward  velocity  results  in  a 
slight  change  of  the  boundary  condition  at  the  hull  surface.  If  we 
consider  the  case  of  head  seas,  the  oscillation  of  the  ship  is  heaving  and 
pitching.  The  effect  of  surging  is  omitted  because  of  the  higher  order. 
Then  the  boundary  condition  on  the  hull  surface  is 


^1  .  ITI,  ,  1T, 

3n’  (Zg-xip-W  nz  ~  U(zg-xi|>)  3n,  (3z  J  ~  9n, 


(276) 


where  n^  =  3z/3n',  n'  being  the  outward  normal  to  the  contour  of  the 
transverse  section.  The  first  term  comes  from  the  relative  velocity  of 
the  section  and  the  last  term  means  the  relative  velocity  of  the  orbital 
motion  of  the  incident  wave.  However,  on  account  of  the  second  term,  the 
boundary  condition  is  not  determined  by  the  relative  velocity  of  the 
section  only.  If  we  assume  U  and  u>  are  both  of  the  order  of  unity,  the 
inner  solution  is  associated  with  the  free  surface  condition  such  as 


=  -  U  C 


w  „  2 
3z 


at  z  =  0 


(277) 


The  source  density  in  this  case  can  be  determined  in  a  manner  similar  to 
the  radiation  problem. 


>(x)  =  -  ^  (iuH-U  {B(x)(zg-xip-Cw)} 


(278) 


In  the  case  of  U  =  0(e  )  and  w  =  0(1),  however,  we  have  to  solve  the 

two-dimensional  problem  with  the  free  surface  condition 


*i.£t  ..vtLh 

3z  g  1  iz2 


(279) 
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Putting 


<t>'  =  4>,  +  U  C 


9^c 

3z 


we  have  the  boundary  condition  at  the  free  surface 


(280) 


=  o 
g 


3z 


(281) 


while  the  boundary  condition  on  the  hull  surface  can  be  written  as 


(zg-xii>-£w)  -  U(zg-x<|j-Sw) 


(282) 


If  we  omit  the  second  term,  we  get  the  usual  form  of  boundary  conditions 
for  a  heaving  cylinder  with  vertical  velocity 

v’(hwh)  <283> 

The  amplitude  of  the  vertical  oscillation  is  then 


Z 


A 


[  V  [  /  a) 


(284) 


If  we  write  the  amplitude  of  waves  generated  by  the  cylinder  as  and 
define  a  two-dimensional  pulsating  source  of  strength,  a  eialt,  which 
generates  the  radiating  wave  of  amplitude  £  ,  there  is  the  relation  for 
the  amplitude  ratio  A  as 


A 


^A  2TTPU) 
ZA  =  *ZA 


(285) 


We  have  the  relation  between  the  density  of  the  three-dimensional  source 
and  that  of  the  two-dimensional  source  as 
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|»  i*m't'**  +*#r' ■•**>■>+■• *"*  ty* 


a  =  2m 


(286) 


Therefore,  we  obtain 


|m(x) | 


=  =  gA[v| 


(287) 


The  determination  of  the  phase  of  m(x)  needs  some  consideration.  Choose 
the  vertical  velocity  of  the  cylinder  expressed  by 


V  =  — —  (A.  cos  ci)t+B_  sin  wt) 
TTOiB  0  0 


(288) 


in  such  a  way  that  the  source  term  has  the  phase  like  sin  a>t.  If  we  write 
the  vertical  movement  of  the  cylinder  as 


Z  =  Z.  e 

A 


itot+ig 


(289) 


the  vertical  velocity  becomes 


io)t+ig+iTT/2 


V  =  a)  Z,  e 
A 


Then  the  phase  angle  of  the  source  is  given  by 


61(x)  =  a  +  g 


where 


tan  a  =  Bq/Aq 


(290) 


(291) 


The  amplitude  ratio  of  the  radiating  wave  is,  on  the  other  hand,  given  by 
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A  =  Z 


2 

A  _  TT  U)  B 

A  2«^R 


(293) 


and  the  absolute  density  of  the  source  is  given  by 


=  gAjvL  . 


(294) 


8/VB  0 


The  coefficients  A^  and  are  determined  by  the  standard  calculation  of  the 
two-dimensional  problem  of  a  heaving  cylinder.  The  density  and  phase  of 
dipoles  which  represent  lateral  oscillations  can  be  determined  in  a  similar 
way.  Since  the  function  U  9<}>q/3z  does  not  include  the  source  term  nor 
horizontal  dipole,  the  density  of  sources  and  dipoles  of  <p'  is  identical 
with  that  of  (J>^.  It  is  not  easy  to  give  a  general  formulation  for  the 
effect  of  the  second  term  of  Equation  (282) ,  but  it  can  be  expressed  by  a 
form  of  a  correction  term  to  the  vertical  velocity  in  the  case  of  semi¬ 
circular  cross  section,  since  the  hull  boundary  condition  in  this  case  is 
given  by  the  form 


=  b&t  (If  +u  It)  {b(x)  : 


(295) 


If  we  assume  the  above  relation  in  arbitrary  cros^  section,  we  can  put  the 
vertical  velocity  of  each  section  in  the  form 


v  ■  TOO  {1M  b) 


(296) 


A  similar  approximation  can  be  applied  to  the  lateral  oscillation. 


Numerical  Examples  for  Added  Resistance 

In  the  previous  sections,  formulae  for  calculating  the  added  resistance 
and  steady  side  drift  force  are  presented.  The  formula,  which  is  originally 
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a  three-dimensional  type,  has  been  simplified  to  a  great  extent  by  assuming 
the  frequency  high  enough.  As  a  numerical  example,  the  model  of  Series  60 
and  C  =  0.7  is  employed.  Here  we  consider  three  kinds  of  calculation 
methods  in  carrying  out  the  numerical  work. 

Method  1:  Assume  U  =  0(1)  and  uj  =  0(1). 

The  added  resistance  is  calculated  by  a  three-dimensional  formula. 

The  boundary  condition  at  the  free  surface  in  near  field  is 

34>, 

—  +  U  X,  — ~  =0  at  z  =  0  (297) 

o  z  w  .  I 

3z 

The  density  of  the  source  distribution  in  head  sea  waves  takes  the  form 

m(x)  =  -  ijjr  ^i(iH-U  (B(x)  (z  -xijj-r^) }  (298) 

Making  use  of  this  source  distribution  in  Equation  (265)  or  (268)  to  calcu¬ 
late  the  Kochin  function,  a  consistent  approximation  in  the  present  case  is 
obtained.  However,  the  result  of  numerical  computation  yields  enormous 
values  of  added  resistance  which  is  hardly  compatible  with  measured 
results.  The  present  formulation  assumes  that  the  disturbance  by  the  ship 
hull  is  represented  by  a  distribution  of  singularities  along  the  axis  of 
the  ship  which  is  taken  on  the  undisturbed  free  surface.  However,  the 
above  mentioned  results  indicate  that  the  singularities  on  the  free  surface 
generate  too  strong  disturbance.  In  order  to  avoid  this  difficulty,  we 
assume  the  singularity  distribution  a  little  below  the  free  surface.  As  a 
mean  depth  of  the  singularities,  let  us  take  the  mean  depth  of  disturbance 
given  by  yT  where  y  is  the  vertical  prismatic  coefficient  and  T  is  the 
draft  of  the  ship.  Then  the  Kochin  function  is  expressed  as 

H*(m)  =  e  m^TJ"m(x)  eimX  dx  (299) 
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1/2 

Method  2:  Assume  U  =  0(e  )  and  w  =  0(1) . 

The  boundary  condition  at  the  free  surface  in  near  field  is 

34>  2  32<j>_ 

3T "  r  *1  +  u  'w  = 0  at  2  =  0  (300> 

The  source  density  is  determined  from  the  two-dimensional  solution  of  a 
heaving  cylinder.  A  consistent  approximation  for  the  added  resistance  is 
obtained  by  the  three-dimensional  formula.  Equation  (266). 

1/2  -1/2 

Method  3:  Assume  U  =  0(1)  or  0(e  )  and  to  =  0(e  ) . 

The  boundary  condition  at  the  free  surface  is  the  same  as  for  Equation 
(300)  in  Method  2.  The  source  density  is  determined  from  the  two- 
dimensional  solution  of  a  heaving  cylinder  with  the  hull  boundary  condition 
of  Equation  (296)  which  includes  the  effect  of  the  forward  speed.  The 
added  resistance  is  calculated  by  the  simplified  formula  of  Equation  (272), 
but  the  solution  is  inconsistent  in  the  order  of  approximation. 

Results  of  computations  of  the  added  resistance  of  the  Series  60 
model  by  means  of  various  methods  at  several  Froude  numbers  are  illustrated 
in  Figures  11  through  13.  Results  of  experiments  are  also  shown  for 
comparison.  It  is  observed  that  the  best  agreement  between  computations 
and  measurements  is  obtained  by  the  calculation  by  means  of  the  simplest 
formula.  Method  3,  while  more  consistent  methods  can  provide  only  less 
accurate  predictions. 

CONCLUDING  REMARKS 

Readers  may  already  have  recognized  that  the  four  topics  discussed  in 
this  report  have  been  arranged  in  the  order  of  simplicity  of  their  physical 
phenomena.  However,  what  is  beyond  our  expectation  for  us  to  find  out  is 
the  fact  that  the  seemingly  simplest  problem  such  as  the  wave  resistance  in 
the  uniform  motion  presents  the  greatest  difficulty  in  attaining  satis¬ 
factory  agreement  between  theoretical  computations  and  measured  results, 
while  theoretical  predictions  for  much  more  complex  phenomena,  such  as  the 


Figure  11  -  Added  Resistance  Coefficient  of  Series  60 

Model  for  C.  -  0.7  and  f  =0.1 
b  n 


Kaw 


X/L 

Figure  12  -  Added  Resistance  Coefficient  of  Series  60 

Model  for  C,  =  0.7  and  f  =0.15 
b  n 
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Figure  13  -  Added  Resistance  Coefficient  of  Series  60 

Model  for  C,  ■  0.7  and  f  *=0.2 
b  n 

added  resistance  in  waves  can  show  much  better  agreement.  Furthermore, 
what  is  really  surprising  is  that  the  most  rational  formulation  brings 
sometimes  most  unrealistic  results,  while  much  simpler  approximations  can 
often  result  in  much  better  predictions.  This  is  the  attribute  of  engi¬ 
neering  and  the  fun  of  engineering.  Mathematics  has  nothing  to  do  with 
this  respect. 

Concluding  this  lecture,  I  wish  to  express  my  sincere  hanks  to  Mr. 
Justin  McCarthy  for  his  kind  arrangements  which  enabled  me  to  have  this 
opportunity,  and  to  Dr.  Kwang  June  Bai  and  Dr.  Ming  Shun  Chang  for  taking 
care  of  me  every  day.  My  thanks  also  to  Dr.  Allen  Powell,  Dr.  William 
Morgan,  Dr.  Wen-Chin  Lin,  Dr.  Bohyun  Yim,  Dr.  Choung  M.  Lee,  Dr.  Thomas 
Huang  and  all  friends  at  DTNSRDC  for  their  warm  hospitality. 
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I  should  be  happy  If  this  lecture  series  can  contribute  to  everlasting 
friendship  between  Japan  and  the  United  States. 
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APPENDIX  A  -  OSCILLATION  OF  A  SHIP  WITH 
SIX  DEGREES  OF  FREEDOM 

Consider  a  ship  floating  on  a  uniform  stream  and  making  oscillations 
of  small  amplitude  about  a  fixed  point.  Take  the  x-axis  in  the  direction 
of  the  uniform  stream  and  the  z-axis  vertically  upwards.  The  fluid  motion 
around  the  ship  is  expressed  by  the  velocity  potential  of  the  form 


<j>  =  Ux  +  U*Q  +  *1 


(301) 


where  U  is  the  velocity  of  the  uniform  flow  and  c()^  gives  the  periodical 
disturbance.  Since  the  amplitude  of  oscillation  is  small,  the  steady 
potential,  U(J>q,  is  determined  from  the  boundary  condition  when  the  hull  is 
fixed  at  an  average  position  and  identical  with  the  velocity  potential  for 
the  uniform  motion  of  the  ship.  Now  we  consider  the  coordinate  system 
Xq,  y^,  Zq,  which  is  fixed  to  the  ship  with  the  origin  at  the  center  of 
gravity.  The  coordinates  of  the  center  of  gravity  at  the  average  position 
are  written  as  xG,  y^,  zG>  and  the  displacement  of  the  center  of  gravity 
as  Designating  the  angle  of  rotation  about  each  axis  by 

0^,  02,  0^,  and  omitting  higher  order  terms,  one  can  write 


x0  =  X  -  XG  "  h  +  03(y-yG)  "  e2(z'ZG) 
y0  =  y  -  yG  ‘  K2  +  ®i(z-zG)  -  e3(*-*G> 
Z0  =  2  "  ZG  "  h  +  02(X-XG)  "  0l(y"yG) 


\ 


> 


/ 


(302) 


If  the  hull  surface  is  expressed  by  the  equation 


F(x0’y0’z0)  =  0 


(303) 


the  boundary  condition  on  the  hull  surface  is  given  by  the  equation 
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(If  +  +x  h  *  *y  fe  +  *2  If )  F(xo'!'o>zo) 

-  -  K1-*,(nc)«!(-c)l  Fx0  -  (E2-91(*-'g)+93<*-*g)>  Fyo 

-  U3-e2o-*G)«i(y-yG))  fZq  ♦  \<FXo-03Fyo+Vzo> 

+  W+Vx.AV  +  *z(Fz0-92Fx0+eiFy0>  "  0  (3M) 


0  0  0 


where  the  dot  means  the  differentiation  with  respect  to  t.  If  we  write  the 
direction  cosines  of  the  outward  normal  as 


„i--p  M,  *2--r  M,  »3  =  -Vi 


(305) 


where 


+  F 


A=  y/V  +  F2 

y  xo  yo 

the  boundary  condition  on  the  hull  surface  becomes 

{C1-03(y-yG)+e2^z'zG^  ni  +  ^2-91(z_zG)+e3(x_xG^  n2 

+  {C3-02<x_xG)+ei^y_yG^  n3  "  ^x(ni-03n2+02n3) 


-  4>y (n2-t-0 3n1-9-Ln3)  -  4>z (n3"02n1+01n2)  55  0  <306) 

Now  we  write  the  steady  portion  of  the  fluid  velocity  as 

u  *  U  +  U94>q/8x,  v  “  U34>q/ 3y »  w  “  Ud<p^/dz  (307) 

and  define  vectors  VQ  -  (U,0,0)  and  V  *>  (u,v,w).  If  we  employ  vector 
notations 


I 


i 

\ 
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•w 


§  -  (CL,C2,C3),  e  =  (01,e2,03) 

- 

r  -  (x-xG,y-yGz-zG) 

the  fluid  velocity  at  the  moving  hull  surface  is  expressed  by 

V  =  y0  +  [(|+0xr)V]  V 


(308) 


(309) 


Substituting  Equation  (301)  for  <J>y,  4>z  in  Equation  (306)  and  making  use 
of  Equation  (309)  together  with  the  boundary  condition  for  4>G  given  by 


n1^0x  +  Voy  +  n3*0z  +  ni  =  0 


(310) 


and  the  irrotationality  condition 


u  =  v  , 
y  x 


v  =  w  , 
z  y’ 


w  *  u 
X  z 


(311) 


we  obtain  eventually 

•  •  •  • 

^lnl  +  ^2n2  +  £3  n3  +  ei  {(yyG)n3-(z-zG)n2} 

+  9 2  {(z-zG)n1-(x-xG)n3>  +  03  { (x-x^nj-ty-y^n^ 

-  E, 1  3u/3n  -  ^2  9v/3n  -  £3  3w/3n 


-  0X  ^  ((y-yG)w-(z-zG)v}  -  02  {(z-zG)u-(x-xG)w} 

3 

-  03  ^  {(x-xG)v-(y-yG)u}  =  (312) 


This  is  the  boundary  condition  for  the  periodical  potential  <f>^.  Then  the 
periodical  potential  is  constituted  by  12  components  such  .is 


♦l 


*  «1*1 
+  Vl 


+  k2*2  +  £3*3  +  +  ^5  +  03*6 

+  E2'P2  +  £3*3  +  e^4  +  02^5  +  03^6 


(313) 
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Boundary  conditions  which  are  satisfied  by  each  component  on  the  hull 
surface  are 


9^/80  =  n^  9^2/9n  =  n2,  9^>3/9n  =  n^ 

9^4/9n  =  (y~yG)  n3  -  (z-zG)  n£  E  n4 

9^5/9n  =  (z-zG)  n1  -  (x-xG)  n3  E  n5 

9v?6/9n  =  (x-Xg)  n2  -  (y-yG)  E  n& 

94<^/3n  "  -  9u/9n  E 

9^2/9n  -  -  9v/9n  E  m2 

9^3/9n  =  -  9w/9n  E  m3 

9i^4/9n  =  -  |^  {(y-yG)w-(z-zG)v}  E  m4 

dip5/dn  =  -  {(z-zG)u-(x-xG)w)  E  m5 

9V9n  =  ~  h  {(x-xG)v-(y-yG)u}  =  m6 


/ 

Forces  and  moments  acting  on  the  hull  are  expressed  by  the  integral 


F 


i 


p  dS 


n  =  1,  2,  3 
n  =  4,  5,  6 


for  force 
for  moment 


(314) 


(315) 


Omitting  terms  of  the  second  order  with  respect  to  the  oscillation  ampli¬ 
tude,  the  fluid  pressure  is  given  by 


where 


P  -  Po 


,2  s*o  .  u2 


pg  is  the  hydrostatic  pressure. 


34) 

9t 


1  3<f>5  3$,  3<t>, 

-  +  u^r  +  v  9r  +  w  ir  (316> 

We  divide  the  above  pressure  into 


■  -  p"2(^r  +  l  lv+0l2) 

/94>  90x 

~ '  p\ar 411 


30,  301 

■+V  -  1 "  4-w  - 

9y  9z 


(317) 


On  Integrating  the  pressure  over  the  moving  hull  surface,  a  periodical 
term  comes  from  the  integration  of  p^  because  of  the  periodical  variation 
of  the  hull  surface.  We  regard  this  as  a  correction  to  the  restoring 
forces  due  to  the  forward  motion.  Forces  and  moments  due  to  the  period¬ 
ical  potential  are  given  by 


'i  ■  °JI 


dS 


(318) 


where  S  is  the  hull  surface  at  the  average  position  and 


01  =  ty^/dt 


APPENDIX  B  -  THEOREM  OF  TUCK  AND  RELATIONS 
DERIVED  THEREFROM 

Tuck  has  proved  a  theorem  with  respect  to  the  integral  in  Equation 
(318)  as  follows 

Tuck's  Theorem 

If  ()>  is  a  harmonic  function  defined  in  the  lower  half  space  outside 
the  surface  SQ  which  has  a  vertical  tangent  at  z  =  0,  and  V  -  (u,v,w)  is 
the  velocity  of  an  irrotational  motion,  which  satisfies  the  boundary 
condition 


un^  +  vn^  +  wn^  *  0 


(319) 


on  SQ,  the  following  relation  is  valid. 


jj  [m^<j>+n^(y*V(}))  ]  dS  =  -  j  n^wds 


(320) 
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where  m^  is  defined  in  Equation  (314)  and  Lq  is  the  intersection  of  Sq 
with  the  plane  z  =  0.  This  theorem  can  be  proved  by  means  of  Stokes' 
theorem  as  Tuck  did,  but  here  let  us  derive  the  same  result  by  integration 
by  parts. 

In  the  first  place,  we  put  i  =  1  and  consider  an  integral 


[i  *JL  <UV 


+v4>  +w<}>  )  ndS 
y  z  1 


-// 


sgnx(u4>x+vi})^+w(j)z)  dydz 


If  we  write  the  equation  of  as 


x  =  f(y,z) 


we  have 


3y  (4>)x=f(y,z)  ^y  +  ^x  3y 


Integrating  by  parts 


JdzJ vtj)ydy  = J dz  J  |v 


9Wx~f  .  3x 
3y  V<Px  3y 


dy 


=  "Jdz/ (V**  If)  -j d2/ V^x  If  dy 

=JJ  (vy<J>n1-vx0n2-v<})xn2)  dS 


By  a  similar  way 


.  3x 

3z  WK  3z 


J*dyJ*  w<J>zdz  =  J*dy J*  jw 

'  Jdy  •  v<Dz=0  -  Jdy  J  (wz+vx  ff  )  *dz  -  Jdy  J w <fx  |f 


dz 


-  J  n^wds  -  JJ(wz<})n1-wx<|)n3-w<f)xn3)  dS 
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Adding  the  above  equations  and  making  use  of  relations  of  continuity  and 
irrotationality, 


V  +  w  =  -  U  ,  V  =  u,  w  =  u 
y  z  x  x  y  x  z 


together  with  Equation  (319),  we  obtain 
I 


1  =ff  (uxI1l+Uyn2+Uzn3)  *dS  '  I  nl4>wds 

Lo 

= rr  ^ds  ~  /  ni^wds 

Lo 

=  -jJ^m^^dS  -J  n^^wds 


A  similar  relation  is  obtained  when  i  =  2, 


X2  =  - 


JJm2<j>dS  -  n. 


(}>wds 


In  the  case  of  i  =  3,  we  write 

I. 


II 


(u<|!x+v<}>y+w<J>z )  dxdy 


We  get  similarly,  as  before 


-jju^dxdy  =  -  JdyJ  (u  f±  -u(j)z  |f) 

=  J/(VUz  10 


dx 


If 


(f>dxdy  +  I  I  u<f>z  dydz 


=  'if  <U^n3"Uznl)  +Uu4>Znl 

-JJv^dxdy  =  -  Jdxj  (v  f*  -v<j>z  fj)  dy 

“JJ  If)  4,dxdy  +U  v^zdxdz 


Summing  up  the  above  equations,  we  obtain 


=  I  I  (w  n.+u  n,+v  n„)  4>dS 

JJ  2  3  2  1  z  2 

-JJ  s  *ds 


m^(})dS 


Then  the  relation  of  Equation  (320)  follows  when  i  =  1,  2,  3.  It  can  be 
shown  that  a  similar  derivation  is  applied  to  the  cases  of  i  =  4,  5,  6, 
if  we  introduce  the  definition  of  n^,  n^,  n^  in  Equation  (314). 


Reverse  Flow  Theorem  and  Haskind's  Relation 

Applying  the  above  theorem  to  Equation  (318)  and  making  use  of  the 
relations 

9$  3i|> 

GhT  =  V  3iT  *  mi 


one  can  express  the  forces  and  moments  in  the  form  like 

io)t 


Fu  ’ e 


5ijf  ‘‘“W  I;  (1"VV  ds 


(321) 


where  the  simple  harmonic  displacement  is  expressed  by  eiaJt  The  line 

integral  term  vanishes  because  we  have  assumed  that  w  =  0  on  the  plane 
z  =  0.  Now  we  ass 
condition  such  as 


z  =  0.  Now  we  assume  that  and  ^  satisfy  the  linearized  free  surface 


a  *  aV  a$ 

-  A +  2iwU  air +  u  r-r  +  8  ar  = 0 

3x 


a^  aV  3$. 

-  «  *±  +  2iu>U  ^  +  U  +  g  _  =  0 

3x 


(322) 


and  the  radiation  condition.  Next  we  define  velocity  potentials  4> 

*  1 
and  which  satisfy  the  same  boundary  condition  on  the  hull  surface  as 


those  of  4^  and  when  the  uniform  flow  has  an  opposite  direction. 
Then  there  are  relations 


34>i  94>. 


9n  9n  ’  9n 

while  the  free  surface  conditions  are 


9^1 

9n 


*  2  *  * 

3$  9$  3<f» 

2  *  1  2  i  i 

-  oo  4>  -  2i(jU  -r —  +  IT  — +  g  =  0 

i  9x  „  2  °  9z 

9x 


*  2  *  * 

2  *  S»!  2  S  S1  3*l 

rr  +  *  sT  ■ 0 

9x 


If  we  put,  for  simplicity. 


*  *  * 

iu)4>^  +  4^  =  4>^,  10)4^  +  = 


(323) 


(324) 


(325) 


the  expression  for  forces  and  moments  takes  the  form 

* 

rr  af 

iiot 


ru  • e 


p£J/s 

b0 


94> 


9n 


dS 


(326) 


We  apply  Green's  theorem  in  the  space  bounded  by  S,  E^,  E,  and  Eg  such  as 
illustrated  in  Figure  3. 


If. 


WMB 


l  3*i  * 

I  — J-  _  - — -  d> 


9n 


$ 

9n  j 


dS  =  0 


If  we  consider  the  case  of  infinite  depth,  the  integral  on  E  vanishes, 

D 

while  the  integral  on  E  vanishes  as  the  radius  of  the  cylinder  tends  to 
infinity.  The  integral  on  the  horizontal  plane  is 

-//  (*1  l/'  ST  *j)d’“ly 

^0 


115 


116 


where  is  the  line  integral  in  Equation  (325) .  If  the  waterline  of  the 
ship  is  slender,  becomes  small.  By  neglecting  we  have  the  reverse 
flow  theorem 


T  .  =  T., 
iJ  Ji 


(331) 


This  reciprocal  relation  was  first  shown  by  Hanaoka  in  the  case  of  a  thin 

ship.  Timman  and  Newman  gave  a  proof  to  the  general  validity  for  non- 

slender  ships,  with  a  somewhat  intuitive  hypothesis. 

When  a  train  of  regular  waves  is  superimposed  on  the  uniform  flow, 

the  periodical  potential  is  composed  of  the  incident  wave  potential  <|> 

and  the  diffraction  potential  if  the  ship  is  fixed  in  the  stream.  The 

wave  excitation  forces  and  moments  are  given  by  Equation  (318)  if  <|>  is 

substituted  by  d)  +  d>„,  or 
w  D 


pJJ  <  W 


(332) 


If  we  apply  Green's  theorem  as  before  and  make  use  of  the  boundary 
condition 


3*w  % 

— —  +  —  =  0 

3n  3n 


(333) 


on  Sq,  we  can  write 


■ '  CJI  (♦. 


S*w  * 
3n 


dS  + 


(334) 


If  the  line  integral  is  omitted,  the  result  is  similar  to  the  relation 
which  was  given  by  Haskind  in  the  case  of  zero  forward  speed.  The 
integral  on  the  large  cylindrical  surface  E  is  related  to  the  Kochin 
function  for  the  radiation  problem.  Therefore,  the  exciting  forces  and 
moments  are  derived  from  the  solution  of  the  radiation  problem  in  still 
water.  It  should  be  noted  that  the  above  formula  does  not  give  a  con¬ 
sistent  expression  for  the  wave  excitation  if  the  finite  forward  speed 
is  present  as  mentioned  in  the  section  on  wave  pressure  on  slender  ships. 
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